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ABSTRACT 

Using  the  Wlgner-Moyal  density  and  operator 
techniques,  W.K.B.  like  solutions  are  obtained  in 
closed  form  for  a  system  of  non-interacting  particles 
in  an  arbitrary  potential  of  arbitrary  dimension. 
For  three  dimensions  the  first  three  terms  of  the 
asymptotic  expansion  give  the  Thomas  Fermi  approxima- 
tion with  the  WeizsScker  correction  reduced  by  l/9  plus 
an  oscillatory  term  of  amplitude  6  times  the  reduced 
WelzsScker  correction.   For  one  dimension  the  first  two 
terms  correspond  to  a  result  of  Kohn  and  Sham. 


Introduction  ' 

In  the  late  twenties  when  the  Thomas -Fermi  model  for 
atom  was  first  Introduced,  It  was  recognized  as  being  crude. 

During  the  early  thirties  various  ad  hoc  corrections 
were  attempted.   Of  these  the  most  prominent  were  the 
Dlrac  exchange  term  due  to  P.  A.  M.  Dlrac  [l]  and  the 
Welzsacker  term  due  to  C.  P.  van  Welzsacker  [2]. 

At  the  same  time  the  Hartree  and  Hartree-Fock 
approximations  emerged.   These  approximations  give  highly 
satisfactory  results,  but  the  computations  are  lengthy  and 
only  feasible  for  atoms.  Ions,  or  molecules  either  with  a 
very  small  number  of  electrons  or  with  a  small  number  of 
electrons  In  a  highly  symmetric  set  of  orbits. 

However  the  success  of  the  Hartree  and  Hartree-Fock 
approximations  led  people  to  seek  some  understanding  of 
the  connection  between  those  approximations  and  that  of 
the  Thomas-Fermi  model. 

In  1957  in  the  work  of  A.  S.  Kompaneets  and 
E.  S.  Pavlovskll  [3],  D.  A.  Klrzhnlts  [4],  and  S.  Golden  [5] 
it  was  shown  that  If  one  expanded  the  diagonal  part  of  the 
Hartree-Fock  density  matrix  In  Increasing  powers  of  -^ 
that  the  zero  order  term  was  the  Thomas-Fermi  result . 
Proceeding  along  this  way  the  same  men  found  that  the 
first  correction  to  that  (the  -fi  term)  gives  the  Welzsacker 
term  reduced  by  I/9  plus  the  Dlrac  exchange  term. 


However  the  resulting  computations  were  unsatisfactory 
and  it  was  observed  that  this  approach  would  never  give 
the  oscillatory  behavior  corresponding  to  the  shell 
structure  of  Hartree  and  Hartree-Fock  solutions. 

More  recently  Von  W.  Macke  and  P.  Rennert  [6]  using 
operator  techniques  and  W.  Kohn  and  L.  J.  Sham  [7]  using 
Green's  functions  have  successfully  remedied  this  last 
problem  for  one  dimension  with  very  encouraging  results. 

In  this  paper  we  will  show  how  formulating  the 
required  equations  in  terms  of  the  Wigner-Moyal  density 
makes  it  possible  to  generalize  the  one  dimensional 
results  of  Von  W.  Macke  and  P.  Rennert  to  D  dimensions. 
Since  our  approximations  break  down  near  r  =  0  for  a 
Coulomb  potential,  we  will  also  discuss  boundary  conditions 
at  that  point. 


Part  I 

Section  1.   Generalized  Remarks. 

In  a  paper  published  In  I963,  W.  Macke  and  P.  Rennert 
[6]    following  a  technique  of  D.  A.  Klrzhnlts  [k]    and 
S.  Golden  [5]  were  able  to  write  down  the  Dirac  density 
In  terms  of  an  operator  integral  for  a  many  body  system 
In  which  all  particles  experience  the  same  potential. 

Prom  this  Integral,  Macke  and  Rennert  showed  that 
the  Thomas-Fermi  solution  dropped  out  as  a  zero  order 
approximation  In  the  D  dimensional  case.   And  then  restrict- 
ing themselves  to  potentials  either  one  dimensional  or 
radially  symmetric,  they  were  able  to  get  W.  K.  B.   like 
solutions  as  first  order  approximations. 

In  this  paper,  we  will  show  how  the  results  of  Macke 
and  Rennert  can  be  generalized  to  the  case  of  a  general 
D  dimensional  potential. 

It  is  hoped  that  our  results  will  make  it  practicable 
to  get  satisfactory  computations  for  binding  energies  of 
heavy  molecules. 

We  now  begin  In  the  same  fashion  as  Macke  and  Rennert. 

For  a  system  of  N  ==  2A  independent  particles  in  a 
common  D  dimensional  potential,  we  have: 

(I. 1.1)         H(r)  ^  (r)   =  E  ^  (r) 

'v>   n  '^'       n  n  'V' 

where 


^2 

The  Dlrac  density  for  the  lowest  energy  state  Is: 

A=N/2        ^ 

(1.1.2)  P(r,r-)   =   2  JZ  V'^C^)  V^^C^')  . 

n=l 

The  factor  2  was  Inserted  above  to  take  care  of  spin. 

Since  the  ii   (r)  form  a  complete  system 
n  ~ 

00  ^ 

(1.1.3)  ZZ:  ^^(r)  ^^(rM   =  5(r-r')  . 

We  can  now  rewrite  (1.1.2)  In  the  form: 


OP      ,  ^ 

(1.1.4)   p(r,r')  =  2^0(A^-n)  ^^(r)  ^^(r') 

n=l 


where 


and 


1   if  a  >  0 

0(a)  =  <  1/2   a  =  0 

0     a  <  0 


A  <  A"*"  <  A+  1 


If  we  assume  p(r,r')  is  composed  of  the  lowest 
energy  levels  and  E^  <  E^^]_  then: 


00  ^ 


(1.1.5)   P(r,r')  =  2^0(E-E^)  ^^(r)  ^^(r') 

n=l 

where  E^  <  E  <  E^_|_^. 

But 


f(H)    Tp    (r)    =    f(E    )    ^    (r) 


so 


(1.1.6)  p{r,r')    =   2  y~  0(E-H)?//    (r)^    (r')    =   20(E-H)5(r-r 

Now  we  may  write: 

+00  +le 

(1.1.7)  ®(«^    =   2?1         J  ^U"  ^"^    ' 

-00  +le 

and 

Ip. (r-r ' ) 

(1.1.8)  5(r-r-)=^^e         ^  dp, 

where  h  =  2Trft. 

Combining  (1.1.6),  (I.l.?)  and  (1.1.8)  we  get: 


+ao+le    .  f-r>  "^y  Ip'(r-r') 

-ia)(E-H)        K,    ^r^  -^    ' 


(1.1.9)   p(r,r')=-^^—    / do)  //  e    '^     dp 

h  2Tri    J  ^  JJ 

-oo+ie 


Section  2.   Zero  Order  Approximation. 

Macke  and  Rennert  pointed  out  that  to  get  the  zero 

order  approximation  we  assume  that: 

_w.2   D    -.2 
(1.2.1)      T  =  -^    2^  -2^  and  V(r)   commute. 


From  (1.1.6)  and  (1.1.8)  we  get: 


(1.2.2)    p(r,r')  =^  JJ    0(E-H(p,?))  e 


Ip- (r-r ' ) 


h 


dp 


Now  using  approximation  (1.2.1)  and  then  setting  r' 
we  get  the  zero  order  approximation  for  the  space 
density  p (r) . 

20->   ^  2    D— 1 

(1.2.3)    P(l)=-^     fe^E  -  V(r)  -  ^)  p    dp 


where 


(1.2.4) 


0 


O,  =^2 


O2  =  2Tr 


O^   =   hir       and    in  general 
n^  ==   27r^/2/r(D/2)     . 


Integrating  (I. 2. 5) we  have; 


/■ 


20^  D/2 

-  [2m(E-V(r))]     if  E  >  V(r)  ; 


(1.2.5)      p(r)  =< 


h^.D 


0 


if  E  <  V(r) 


L 


For  the  case  D  =  3,  (I. 2. 5)  is  precisely  the 
Thomas  Fermi  or  T.F.  approximation.   One  only  needs  to 
add  the  requirement  that  V^V  =  -hire   p(r)   to  get  the  usual 
equation  for  the  potential  of  the  T.F.  atom  [8]. 


Section  3.   The  N  Particle  Wlgner-'Moyal  Density 
for  D  Dimensions . 

We  shall  begin  by  rewriting  equation  (1.2.2) 

Ip • (r-r ' ) 
(1.2.2)   p(r,r«)  =  -^  IT  Q{E-E{p,r))    e         ^     dp 

but 

Ip • (r-r ' ) 


O.*    <^  '\^ 


(1.3.1)  e(E-H)e    ^ 


^  (p.+p.)2  ^  -^  ^ 


/T  -1^   ^  p .  +p .  j 

i/  d^"5(r:-r")  Q(E-  YZ        \^^        -   V(r"))  e 


tr 

.  5;;^ -  v^r  ;  j  e 

1=1 

where : 


^  ip-r/h 

p^  acts  only  on  f(r)   (In  this  case   e  '^  "^   ), 

p.  acts  only  on  f(r")  (in  this  case  V(r"),  and 

p.  Is  the  1-th  component  of  p. 

Now  It  Is  clear  we  may  replace  the  operator  p. 
in  (1.3.1)  by  the  corresponding  scalar  p.,  giving  us 
the  equation: 

lp-(r-r  ' ) 


(1.3.2)  0(E-H)e    ^ 


^"    .     ^  A^     VO   'V 


IT  D  (p.+p  ^'^ 

=  //  dr"5(r-r")  0(E-  y-    '^^   ^±'      -   v(r"))  e 


ip- (r-r ' ) 


-h 


or  after  integrating: 


Ip'(r-r')    Ip'(r-r') 
(1.3.5)  ©(E-H)  e  =e         ^  Q(E-  y~       t      ^     -V(r))-1. 

Thus  (1.2.2)  may  be  rewritten  in  the  form 

(I.3.M  p(r,r.)  ^  //  dp  e    ^      -^  0(E-  ^  -^^   -V(r)M. 

Now  from  equation  (A. 1.34)  In  Appendix  1  It  Is  clear 
we  now  have  a  formula  for  the  complex  conjugate  of  the 
Klrkwood  density, 

(1. 5.5)      PictP<£'  °  ^e(E-Z:   V   -  V(r))-1 

Using  equation  (A.I.38)  and  keeping  in  mind  that 
the  Wigner-Moyal  density  is  real  we  get  a  formula  for 
Wlgner-Moyal  density. 


(1. 3. 6)      p(p,r)  =   exp  (^  |3  ^^tI^)  Pk^P'^^ 


and  now  from  (1. 3. 5)  noting  that  j  -5 —  =  p.  , 


If  we  were  to  make  the  T.F.  approximation  that 


8 


dV{r)/dx.  ^  0  I.e.  V(r)  ^  a  constant,  then  (I. 5. 7)  becomes 


(I. 3. 8) 


p(p,r)  ^  ^  e(E  -  H(p,r)) 


Section  4.   Expansion  of  Wigner-Moyal  Density  in  Powers  of  fi. 


One  way  to  extract  something  useful  out  of  (I. 3. 7) 
Is  to  expand  p(p,r)  In  Increasing  powers  of  h,  keeping 


A    -fi  d 


In  mind  that  Pj_  =  i  (j^ 

1 
Doing  this  we  get: 


,D 
(1.4.1)  %r   p(P,r)  =  0(E-  H) 


+ 


25m 


D 


1  1>  J     -'-   J 


1=1 


w 


-  3  IZ  (^)  e^^^E-  H)]  +  O(fi^)  , 

1^    ax^ 

D_  pf 
here  H  =  XI  p^  +  V(r)  . 


1=1 


If  we  Ignore  the  second  order  derivatives  and 
Integrate  out  momentum  we  get  the  space  density 

^    -  ,D/2 


(1.4.2)   p(r)  =<^ 


20p[2m(E-V(r))]  ^   .  ,  D-fD-2)(D-4)    ^^^\T^\^     -, 

^^      [2m(E-V(r))p 


h^D 


0 


E  >  V 
E  <  V 


In  1957^  Kompaneets  and  Pavlovskll  [3]  were  the 
first  to  compute  a  similar  series  for  p(r)  In  powers 
of 'fi.   However  their  method  Involved  using  a  system  of 
hierarchy  equations  with  non-unique  solutions. 

Shortly  afterward  In  the  same  year  D.  A.  Klrzhnlts  [4] 
was  able  to  apply  some  operator  techniques  he  developed  to 
put  such  a  series  for  p(r)  on  a  sounder  footing. 

Although  D.  A.  Klrzhnlts  wrote  down  the  first  few 
terms  for  p(r)  and  not  for  the  more  useful  Wlgner-Moyal 
density  p(p,r),  he  was  still  able  to  use  his  methods  to 
compute  any  expectation  value.   At  any  rate  It  Is  only 
a  short  step  from  his  work  to  equation  (1.4.1)  and  we 
have  used  methods  very  similar  to  his  to  get  It. 

In  i960,  G.  A.  Baraff  and  S.  Borowltz  [9]  developed 
a  very  clean  method  for  getting  a  similar  expansion  for 
the  Green's  function  corresponding  to  the  Wlgner-Moyal 
density. 

The  first  term  In  (1.4.1)  Is  of  course  the 

2 

Thomas  Fermi  approximation.   The  first  term  of  order -tl 

in  the  square  brackets  is  the  source  of  the  second  term 

In  (1.4.2)  which  is  generally  referred  to  as  the 

1/ 
Weizsacker  correction. 

Kompaneets  and  Pavlovskll   pointed  out  that  the 

2 
Dirac- exchange  term  is  also  of  order  -h  and  should  be 

—  Our  "Weizsacker"  term  is  that  of  Kompaneets  and 

Pavlovskll  i.e.  I/9  of  the  original  one. 
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treated  as  a  correction  to  the  T.P.  approximation  of 
the  same  order  as  the  Weizs^cker  term. 

One  could  continue  expanding  p(p,r)  out  further 
in  this  fashion.   However  all  correction  terms  are 
functions  of  the  p. 's  and  derivatives  of  V(r)  multiplied 
by  higher  order  derivatives  of  O(E-H)  which  are  in  turn 
derivatives  of  5(E-H) .   So  it  is  clear  that  additional 
terms  may  possibly  improve  some  expectation  values  but 
they  certainly  will  not  affect  the  actual  Wigner-Moyal 
density  away  from  the  point  where  H  =  E.   Thus  it  is 
clear  that  not  all  expectation  values  would  be  Improved. 

Having  discussed  the  power  series  approach,  let  us 
return  to  (I. 3. 7)  and  try  a  hopefully  more  fruitful  method 


Section  5.   W.K.B.-llke  Approximation:   ^^V/Sx.Sx.  =  0. 


1      3 


From   (1.1.7)    we  may  write 


D      (p.+p.)2 
(1.5.1)    0(E-  YZ-^^ V(r)) 


i^l  2m 


—  Pi       .      ^  P?+2p^p^ 


1      r     dcD  ^-^  ^--^ 


27ri       ,  /  (D 


e  e 


-c»+ie 
Now   let  us   consider   the   operator 


11 


(1.5.2) 


A+B 


±0^(1 


>\2  A 

D_  p.+2p   p 

^        ^    ^   +V(r)) 


1=1 


2m 


where 


{I.5.3K 


A  =   IcD  V(r) 


D_     p^+2p.p. 

B   =    10.  ZZ    (    ^   o^^   ^) 

1=1 


2m 


-  ^  z:  (^^  -4  +2if^p    ^ 


ICO        V fj^'^         <~>  ,rsAx:  "      \ 

2   +21t^Pi  ^3^) 


2m 


1=1 


^x^ 


Sv      a 


D 


^V 


[B,A]    =  i  [«2(  ^  l!v  ^  2  ^  u^  ^)^  21^  ^  p^   u^3 

1=1    ox.  1=1         1         1  1=1  1 


2m 


(I.5.^X 


1=1     ^^i 


m 


[[[B,A],A],A3    =   0 


Now  we  make   the  approximation   that  B,    [B,A] ,    and 
[[B,A],A]    commute  with  one  another.      Doing  this  and 
using   equation    (A. 2. 12),    we  get 


(1.5.5)      e 


A+B  ^         (B+-|[B,A]+-^[[B,A],A]) 


=      e     •  e 


This  approximation  corresponds  to  assuming  S  V/Sx.^x.=0, 
We  should  note  that  In  comparing  this  approach  to  that 

used  in  the  -last  section,  we  are  in  effect  Ignoring  the  second 
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2 
and  third  terms  of  order  'R   in  the  square  brackets 

of  equation  (1.4.1). 

However  these  terms  are  usually  Ignored  anyway 
and  we  are  no  longer  making  the  approximation  that  ^ 
is  near  zero. 

Collecting  (1.5.2),  (1.5.4)  and  (I. 5. 5)  together 

we  have:        ^^ 

D  p.+2p.p. 


i^(^  "  %^"  "  +V(r)) 


(1.5.6)  e 


1=1    2m 


U.V(r).  i|5«  2=  P.  ^  .  ^  ^  (1^)2] 
'~'    2m  ^— ^  ^1   ax.      6m   4— t-  ^dxT'^  ■' 

1=1       1  x=l    1 


=  e 


From  (1.3.7),  (1.5-1),  and  (I. 5. 6)  we  get 
(1.5.7)   p(p,r)  = 


+°o+ie  -^ZIPi-5^ 

2    1       r  do)    -i^E         1=1   ^1 

e     e 


^D  2Tri   J   00 
-oo+ie 

2m  ^7—Y   ^1  dx.     6m  <—rr    Mx.  '  ' 
1=1      1         1=1    1 


Keeping  in  mind  the  approximation  that  ^^y/^^.^x.  =  0, 

J.     J 

it  can  be  shown  by  induction  on  N  that 
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(1.5.8)      (_I^p^^,     exp   (l.H(p,r)   .  i|^  |Z  Pi  II;; 
^  loD^fi^  ^:^   ,av   ^2^ 

[N/2]      „,     ,    ,,N        N-2k      k  •    2„      D  ;.,, 

=     2 u  V      exp(icDH+  -5—  ;^  p     -^-r- 

k=0      (N-2k) I    kl  ^     1=1      1   °^i 

■      lo^^^     ^;:^    aSV   ,2, 
1=1  1 


where    [N/2]    is   the  greatest   Integer   <_  n/2,   and 

u  = 
(1.5.9) 


2nr  5^  Pi  "53r  +  ~5i       5^  ^^' 


±=i   -^  ^^1      ""'"       1=1   ^^1 


.    3^2     JD        >„   2 

1=1        1 


Noting   that 


^  -  a  ,  1^.     ,  S  M 


1  x -^        d  /I  A        6    % 

(1. 5. 10)  e  "--^  ^     =     Z~" 


N=o  n: 


and 


^1        rgZ2]      NI     (-1)^  N-2kk 

(1.5.11)     Z_  M^       ^ —-^ —     u         V 

N=:0  '■^'        k=0      (N-2k)  :    kJ 


GO            k       00  /    n  \N  N-2k 

V      ^: —       (-1) 


k=0      ^-      N=2k      ^^^    ^^'  ' 

v-u                       r         l-Hco^     ^  „      SV  Icu^-h^  ^  /5V    ^2•, 

'             =      exp    [    -  -g^     ZZ  Pi  ^5r  -  -8i?I-  f^  ^^^    J 

1=1               1  1=1           1 


1'^ 


We   conclude   that 


1     T-^^      ^ 

-2  5^  Pi 


e  (expdcDH  +  -^  ;^  p.   ^  +  -^^  S^^' 


(1.5-12) 


.    5^2        D        N,,      2 

=      exp    (ia)H  +-2l|jr     ?^   ^^5^^    ^ 

1=1  1 


Thus    combining    (I. 5. 7)    and    (I. 5. 12)    we  have 


-00+le 


where  |vv|2  =  ^Z   (U")^- 
1=1  "^^i 
From  (1.5.13)  we  may  Immediately  compute  forms 

for  the  space  density: 

00 

p  D-1 

(1.5.1^)  p(r)  =   O.   /  p(p,r)  p    dp 

0 
and  the  kinetic  energy  density: 

00        2 

(1. 5. 15)  t(r)  =  Op  |p(p,r)|^p^-l  dp  . 

0 

To  get  p(r)  and  t(r)  from  (I.5.13)  we  consider  the  integral; 

°^   J^  N-1        C^      p^  N/2   ^   4^ 

(1. 5. 16)  Oj,  /  e  2m  p    dp  =  ^  (f )    r(|)  e^. 

0 
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Combining    (I. 5. 13),    (1.5.1^)   and    (I.5.16)   we  get; 


IttD 


(1.5.17)    p(r) 


-O^  D/2 

^:(2m)         r(5)    ^ 


+oo+ie 


h 


D 


2'      2Tri 


•CD  +ie 


doo 
D+2 


loD^i 


•la)(E-V(r))+^t|^    |VV 


and   similarly   combining    (I.5.I5),    (I. 5-15)    and    (I.5.16) 
we  get 


(I. 5. 18)      t(r)    = 


-^ 


D 


2m  h 


-1(J0 


T^^o  iTr(D+2)    4<3D+i£ 

,„    V    2     p/D+2x    e r         du3 

(2m)  r(-5-)   — 7^^ /        -^pi 


27ri 


-00+ie 


O) 


T 


(E-V(r)).^    |W|2) 


Equations    (I. 5.13),    (I. 5.1?)    and    (I.5.I8)   may  be  written 
in   somewhat  neater   form. 
From   (1.5.13): 

+oo+le 

do)  e 
-00+ie 


^p(p,r)  p  n 

1.5.19)    5e "  ^  2i 


.    3^2 
-ia)(E-H(p,r))+^flj^    IVV|2 


so   that 
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ap(p,r)     p  -2m(E-H(p,r))  ^ 


where  A. (x)  Is  the  Airy  function.   Assuming  p(p,r) 
then 


=  0 
E=-oo 


^  -2m(u-H(p,r)) 


(1.5.21)  p(p,r)  =  ^   /  7^^-— 27^  A  ( ^ ^>._^)  du 


-00   '    '~ 


2m(u-H(p,r)) 

or,  substituting  t  =  ^ o^r-z.      >    we  get: 

(mfilvvl)'^/^ 


2m(E-H(p,r)) 

^73" 


|mfi|VV|  )■ 
(1.5.22)  p(p,r)   =  ^  y      A^(-t)  dt 


-00 


To  get  somewhat  similar  forms  for  p(r)  and  t(r)  we 
should  first  recognize  that  Integrating  (1.5.14)  and 
(1.5.15)  by  parts  we  have: 


00 


(1. 5. 25)     p(r)  =     -^    J     P^'-S^dp 

0 


and 


00 


0 


27 

—    Handbook  of  Math.  Functions,  edited  by  Milton  Abramowltz 
and  Irene  A.  Stegun  (New  York,  Dover  Publ.,  Inc,  I965), 
p.  447  (10.4.52) . 
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From    (1.5.15) 


^^-^  _l,(E-H(p,r))+if|!|vv|2 


Sp(p,r)  _p    n     T  r 

ti.5.25)    ^^'^^h    J     ^' 

-00+le 
iso   that 

(1.5.26)    -5j^    =    ?"(„«|vv|)2^5    ^     (^Ivvl)'-^/? 

2m[E-V(r)-  |^] 
Substituting   t  =  o7^ —   ,    we  get 


Sp(p,r)  p 

h^ 


(I.  5. 27)  §^      dp  =  4  ^i^-^^       ^^ 


^/^  2m(E-V(r))  wp 

and  p  =    (mR|vv|)^/5( ,       ,"g/.   -   t)!/^ 

(in«  VV   )^/5 

Thus  (1.5.23)  becomes 


2  0,         D/5   r     ^/2 
(1.5.28)   p(r)   =  -^   (mftlvvl)      /  (a-t)    A  (-t)  dt 


-00 


and  (1.5.2^1-)  becomes 


^+2   ^     D+2 


(1.5.29)   t(r)  =  -j^ (mh|VV|)      r(a-t)  ^  A  (-t)  dt 

h^(2m)(D+2)     -        J  1 


-00 
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■where 


(I. 5. 50) 


a 


2m(E-V(r)) 
(nrTilVVD^/^ 


When  D  Is  odd  It  Is  possible  to  integrate  (I. 5. 28) 
and  (1.5.30)  explicitly. 

Prom  equation  (A. 3. 12),  we  see  that: 


D+1 


(-1)   %  a,       ^/, 

(1. 5. 31)      P(r)    -  5 ^   (2iTrtllvV|)^/5 


D  h 


$^V^-7)Vi;!75-i' 


for  D  =  1, 5,5, etc . , 


u=0 


and  from  equation  (A. 5. 11),  we  see  that: 


(1.5.52)      t(r)    = 
D+5 


D+5  D+2 

(-1)    ^   H  Op(2mti|VV|  )    ^ 


d 

du 


(D+2)    2m  h 


D 


for  D  =   1,5,5, 


u=0 


For  the  complete  asymptotic  expansions  of  p(r)  and  t(r) 
for  arbitrary  D  the  reader  should  refer  to  Appendix  ^, 
particularly  (A.i^.20),  (A.J+.21),  (A.H.22),  (A.1+.26),  (A. 4. 2?) 
and  (A. 4. 28). 
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However  in  the  particular  case  where  D  =  3,  with 

2m(E-V(r)) 
a   = >   0 

(mh|VV|)^/^ 


we  find  that 


^   [2m(E-V(r))]    [.-    ^-   ^  ....  - 
5h5         ^  Sa^  ka^  ^ 


(1.5.33)   p(r)  -  -^  [2m(E-V(r))]    [l-  -^-  -2-  sin  |  a^^^] 


(1.3.3^)   tCr)-.  Jt2I_  [2m(E-V(r))]    [  1+  -^ +_1|  cos  |  a^^^] 


and,  with  a  <  0,    we  have 

5/2   ^   5 
(1.5.35)   p(r)  -  -^  [2m(V(r)-E)] 


'^ 


h^       -"  (-a) 


5/.    -^'-"'"' 

(1.5.36)   t(r)  -  -^  [2in(V(r)-E)]    -^ ^ 

4mh5       -"  ^   ^9/2 

(-a) 


These  may  of  course  be  regarded  as  connection  formulas 
Of  particular  Interest  here  Is  the  comparison  of 
(1.5.33)  and  (1.5.35)  with  equation  (1.4.2).   The  first 
two  terms  of  (1.5.33)  are  the  T.P.  approximation  with 
Welzsacker  correction  as  In  (1.4.2).   However,  the  third 
term  In  (I. 5. 34)  Is  an  oscillatory  term  with  an  amplitude 
6  times  the  magnitude  of  the  reduced  WelzsMcker  term. 
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It  should  be  noted  that  this  oscillatory  term  would 
never  have  been  produced  by  the  method  of  expansion 
In  powers  of  -h.   Also  the  method  of  expansion  in  powers 
of  ^  will  never  give  anything  other  than  zero  for  the 
region  V(r)  >  E. 

We  should  also  consider  the  one  dimensional  case 
and  then  compare  our  result  with  that  of  Kohn  and  Sham  [7] 
and  that  of  Macke  and  Rennert  [6]. 

For  one  dimension,  equations  (A. 4. 20)  and  (A. 4. 21) 
give  us: 

J,            1/2     cos  I  p 
(1.5.37)     p(r)  'v.  ^  [2m(E-V(r))]    [l-  ^]  . 

Now  if  we  continue  to  make  the  approximation  that 
V(r)  =  a  +  br  (b  >  0),  then 

_       ^   [2m(E-V(r))]        ^  P    r  ^/^ 

-I  p  =  4  =  Ij   [2m(E-V(u)))    du. 


5  ^     5 


mfi  VV 


fVr 


where  V(rQ)  =  E. 

Also: 

1/2      ^0 


[2m(E-V(r))l  '   ^   1  T  [2m(E-V(u)  )  ]  "'''^du. 


mh  VV 


'\* 


Now,  if  we  substitute  the  parameters  used  by 
Kohn  and  Sham, 
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1/2 
p(E,r)  -  [2m(E-V(r))] 


"0  -1/2 


T(E,r)  =     ^    f     [2in(E-V(r))] 


du 


r  p 


-,   o  0  1/2 

0(E,r)  =  ^      /   [2m(E-V(r))]     du 


then  we  get 

(1.5.58)     p(r).4tp(E,r)  -^ffljflf^] 

which  aside  from  our  spin  factor  of  2  is  precisely 
the  result  of  Kohn  and  Sham. 

Note!   Although  equation  (I. 5. 37)  agrees  with 
the  result  of  Kohn  and  Sham,  it  does  not  agree  with 
the  result  of  Macke  and  Rennert .   It  appears  that 
P.  Rennert  (it  was  his  dissertation)  made  some  error 
in  computing  his  asymptotic  expansion.   In  both  the 
Macke-Rennert  paper  and  Kohn-Sham  paper  graphical 
comparisons  were  made  of  exact  and  approximate  solutions 
for  harmonic  oscillator  and  Coulomb  potentials.   These 
graphical  comparisons  are  a  little  disappointing  in  the 
Macke-Rennert  oaper  whereas  they  were  excellent  in  the 
Kohn-Sham  paper.   The  only  significant  error  on  the 
Kohn-Sham  graphs  appears  for  the  Coulomb  potential 
when  r  is  near  zero.   And  perhaps  this  region  should 
be  attacked  differently. 

■2/   W.  Kohn  and  L.  J.  Sham,  Phys.  Rev.  I37  A  I7OI  (I965) 
equations  (^-l^),  (3.15)  and  (3.16) 

22 


Section  6.   Relation  to  W.K.B.  Method. 

First  of  course  we  have  made  the  approximation  that 

o 

S  Y/hx.bx.   =  0,  which  Is  done  in  the  transition  region 
-^   J 

for  the  W.K.B.  method. 

Now  if  V  =  a  •  r  +  b  where  a  is  a  constant  vector 
and  b  is  a  constant  scalar  and  E  >   V(r)  then 

P   ^/P     p  [2m(E-V(r))]^/2   P  l/2  (-a) 

(1.6.1)  V  I  a^/2^v|  -^ =  1  [2m(E-V(r))]    — ^ 

^   ^  ^  ^  mh  VV         ^         ^  a 


Thus 


3/2    ^   p^  1/2  (-a) 


(1.6.2)  I  a    =  I  J   [2m(E-V(u))]    -^-^  •  du   , 

where  r^  is  a  value  such  that  V(r^)  =  E. 

For  an  attractive  potential  in  one  dimension, 
+a  =  I  a  I  and 

P  3/2     or  ^/2 

(1.6.3)  I  a     =  I  j    [2m(E-V(u))]     du  . 


For  E  <  V,  and  proceeding  similarly,  we  get: 

r 

P     3/2     p   p'  1/2  a 

(1.6.4)  ^  (-a)     =  I   /   [2m(V(u)-E)]    ^^   .  du 

And  here  again  for  a  one  dimensional  attractive 
potential  we  have  the  form 
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3/2    ^   /?  1/2 


(1.6.5)   I  (-a)     ^     i    J      [2m(V(u)-E)]     du  . 


^0 


Of  course  when  VV  Is  not  a  constant  vector  (1.6.2) 
and  (1.6.^)  are  no  longer  Identities. 

If  for  example  we  assume  that  the  potential 
experienced  by  electrons  of  an  atom  at  long  range  is 
-e  /r,    then: 

^T  C   C\  2,   ^''/'^    2   [2m(-E)1     r"^ 

(1.6.6)      ^(-a)    ~  ■% p 

-^  -^  m  ^  e 


whereas 


1/2 
^  1/2       2[2m(-E)]     r 


(1.6.7)      I  f    [2m(V(u)-E)]     du  -- 


^0 


Looking  at  (I. 5. 35)  it  appears  that  for  long  range 
it  would  be  more  appropriate  to  use  approximation  (1.6-7) 
rather  than  (1.6.6). 


24 


Section  7-   Region  Near  Nucleus  of  an  Atom. 

It  appears  that  our  approximation  should  be  good 
In  the  transition  region  where  E  -  V(r)  Is  dominated 
by  Its  linear  part.   Our  approximation  should  also  be 

good  In  the  oscillatory  and  exponentially  decreasing 

2       ^/2 
regions  particularly  If  we  replace  -^  ( +  a. )     by 

the  appropriate  W.K.B.  Integral. 

However  a  glance  at  equation  (I. 5. 28)  shows  that 

o 
we  have  a  disaster  for  r  ~  0.   If  V  =  -  Ze  /r   then 

|VV|  =  Ze^/r^  so  p(r)  =  0(l/r^).   This  singularity  Is 

such  that  the  Integral  for  the  total  potential  energy: 

p(r)  r  dr     would  not  exist. 

0 

One  way  to  avoid  this  difficulty  Is  to  rewrite 
(1.5.28)  for  D  -  3  In  the  form 

R^  5/2  ?:  .    3/2 

(1.7.1)    p(r)  =  ^  [2m(E-V(r))]     /   (1-  ;J)    A(-t)dt 

^  -00 

2   ^/2 
and  replace  -^  a  which  equals  zero  at  r  =  0  by 

pZe^/{-E)  7  2    / 

^   /        [2m(E+^Y~^l    ^^'  which  Is  not  zero  at  r  =  0. 

However  It  should  be  noted  that  the  justifications 

for  a  W.K.B.  approximation  break  down  in  this  region. 

And  Klrzhnlts   for  one  states  that  a  large  amount  of 
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the  error  In  T.F.  total  energy  computations  comes  from 
Just  this  region. 

So  now  let  us  consider  some  other  way  of  handling 
the  situation. 

A 


(1.7.2)  p(r)   =  2  ^H  ^  (r)  ^^(r) 

k=l   '^ 

-Zr/na      .      V 

(1.7.3)  ^^(r)  -     e      °.  0(r  ) 


-%/ 


where  £   Is  the  azlmuthal  number,   n  the  total  quantum 
number  and  a^  =  the  first  hydrogen  Bohr  orbit  =  f\   /me  , 
Thus  we  see  that  near  the  nucleus  p(r)  Is  dominated  by 
the  s  states. 

For  an  s  state: 


*.2   n2         ^2  h±{r)        r.2 
dr  or 


Multiplying  (1.7-4)  by  r  and  then  letting  r  approach 
zero,  the  Important  terms  of  (1.7.4)  become: 

(1. 7. 5)       -S-— fer-  -  Ze^  i^Jr)   =   0  . 


m 


Multiplying  (1.7-5)  by  i^y.i'^)    we  get: 


2 
(1.7.6)    -  ^  ^*(r)  1^  ^^(r)  -  Ze^  ^*(r)  ^^^(r)  =  0. 


Adding  (1.7.6)  to  Its  complex  conjugate  we  have: 

— /   Leonard  I.  Schlff,  Quantum  Mechanics  (New  York, 
McGraw-Hill  Book  Co.,  In'c^^   1955)  P-  ^5- 
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-h^  d 


2   * 


(1.7.7)    -  ^~:   (^i,(r)^(r))  -  2  Ze"  ^j^(r)  ^^(r)  =  0 


Summing  over  all  s  states  we  have 


fi^  h 


(1.7.8)    -^ 


m 


^ 


> 


s  sta 


^  (r)^ir) 
tes 


-  2  Ze^    ZZ    ^^(r)^r)   =  0  , 


s  states 


or 


(1. 7. 9) 


^  _  J: 


s  states 


\ir)    V'^(r) 


I 


s  states 


^l^(r)  ^(r) 


2Z 
-  —  for  r  near  0 , 

^0 


Now  referring  to  (I. 7. 2)  and  (I. 7. 3)  it  is  clear  that 


(1.7.10)   -^p(r)   = 


dr 


-T-    YH       ^v(r)  ^i^(r)   for  r  near  0; 


dr 


s  states 


and 


p(r)  =     >     ^i,(r)  ^^(r)      for  r  near  0 


s  states 


Thus  the  appropriate  condition  to  use  on  p(r) 
for  r  near  zero  is 


where  a^  = 


« Vme^ . 
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Now  let  us  turn  to  the  kinetic  energy  density, 


(1.7.12) 


tL^)  =  //  P^P'£)  fe  ^P 


2m 


On  the  other  hand  taking  the  Fourier  transform 
of  equation  (A. 1.12)  we  have 


iT-p 


(1.7.13)    //  p(P,r)  e     dp 


2X1  K(r+1  i)k(t-  S  j) 


'  n^  /±      r\j 


k^ 


Prom  (1.7.12)  for  D  =  3 


(1.7.1^)    t(r) 


J- 7 
2m  ^ 


1-   a^ 


J=l  SxJ 


IT'P 


p(p,r)  e 


dp 

'V. 


T=0 


Thus  from  (I.7.I3)  and  (1.7-1^) 


(1. 7. 15)   t(r)  =  -  2(lj^)  ^^  r(v2^^(r))^;(r)  - 


-  2y^i.(r)-yK(r)+  ^v(r.)v2^r(r)] 


'V>  k  "^   'v  k  '^^ 


'k' 


but  since 


we  may  write 


-Ii^'^k=  (^-^te^^^k  ' 


(1. 7. 16)   t(r)  =YZ   (Vk^^^^k^^^-^f^^^k^^^^k^^^+li^^k-^^k^ 


k 
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Substituting  V(r)  =  -  Ze  /r,  multiplying  (I.7.16) 
by  r  and  then  letting  r  approach  zero  we  have 


Section  8.   Bohr  Sommerfeld  Rule  for  Hydrogen  Atom. 

If  we  use  V(r)  =  -  ZeVr  then  p(r)  =  2  ^  ^  (r)^^(r) 
a  sum  over  the  first  N  solutions  of  the  hydrogen  atom. 

Also  E,the  T.F.  cut-off  energy  or  chemical  potential 
if  you  wish,  has  the  property  that 


En  ^  E  -=  ^4-1 

Now  from  (I. 5. 33)  and  (1.5.3^) 


r  t  ( r ) 
(1.8.1)    .rJ:     ~l  (Ze^+Er) 


[1  + 


Jlrl 


8a 


3 


+ 


15     2   3/2t 


8a 


-,  cos  ^  a' 
3     3 


[l--^-^sin|a5/2^ 


8a- 


4a 


3 


but 


(1.7-17) 


r  t(r) 


ize2 


r=0 


Equations  (1.7.17)  and  (I. 8.1)  seem  to  suggest 

2  3/2 
that  at  r  =  0  the  W.K.B.  integral  corresponding  to  ^  a 

should  be  an  angle  in  the  second  quadrant.   Thus 
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Ze^A-E) 
o         r  7^2   1/2 

(1.8.2)  ^   /   [2m((+E)  +^)]     dr  =  (2n+A)7r 

0 
where  n  Is  an  Integer. 

Carrying  out  the  Integration  we  get: 

2  2 

(1.8.3)  E  = ^  ^  .  o       where  ^  <  A  <  1. 

2aQ(n+|)^  ^ 

Of  course  to  get  E  more  accurately  we  should 
adjust  E  so  It  satisfies  the  equation: 

/  p(r)  dr  =  N    where  N  Is  the  number  of  electrons 
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Part  II 

Section  1.   Exact  Wlp:ner-Moyal  Density  for  Non-Interacting 
Harmonic  Oscillator. 

When  discussing  approximation  schemes  it  is  nice 
to  have  some  exact  solutions  for  points  of  reference. 
First  they  provide  a  test  of  sorts  for  the  accuracy  of 
any  given  approximation  scheme.   And  secondly  they  may 
provide  starting  points  for  other  approximation  schemes. 
Thus  in  Sections  1  and  2  of  Part  II  we  will  discuss 
solutions  to  a  system  of  N  non-interacting  bodies  in 
a  D  dimensional,  spherically  symmetric,  harmonic 
oscillator  potential. 

In  Sections  3  to  5  of  Part  II  we  will  discuss 
a  system  of  Interacting  harmonic  oscillators. 

Consider  first  the  one  dimensional  harmonic 
oscillator 


2   2 
(II. 1.1)     -|--^+ikx^u  =  Eu  . 

^^   dx 


This  equation  has  normalized  solutions 
(II. 1.2)     u^(x)   =  N^  H^  (ax  )  e  ^ 


^^       Schlff,  Quantum  Mechanics,  p.  64. 
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with   corresponding   eigenvalue  E     =    (n+  -p)ha3c,   cue   =  v/k/m. 


l/2„n^,,l/2  „_,      4        .„,_,, 2 


where     N^  =    (q/tt  '^    2  n  I )    ^      and  a'"   =  mk/h 


Now  since  the  equation; 


(II. 1.3) 


2i4^T^  ^  +2  ^1^1  ^ 
1=1  ox. 

1 


Eu 


is  separable,  it  is  clear  that  (II. 1. 3)  has  normalized 
solutions  of  the  form: 


(II. 1.4)   u 


m^^mg. 


,m^   - 


(x)  = 


D 

TT  N   H   (ax.) 
.   -,  m .  m .    ,] 

D  =  1   J   J     . 


1   2|  |2 


with  eigenvalue 


D 


D, 


E„  =   (  2_  m,  +  ^)  1113^ 


m 


Now  if  we  introduce  the  notation  that 


(II. 1.5)   H^(ax,ay) 


H  (ax)H  (ay)     l/2   p 

,n;;^     =^    (NX(-)H^(ay)) 


and 


(II. 1.6)   fr(ax,ay)  =      YH  \   (ax^,ay^)H^  (ax^jayg) 


~    m-, +m2+.  .  mp.=m  "1 


Then  we  may  write  the  Dirac  density  for  the  complete  set 
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of  states  with  energy  eigenvalue  E  =  (n  +-p)ha3  as 

^2  D/2  -  |!(|x|2+|yl2) 

(II. 1.7)    P^(x,y)  =  (^)    H^(ax,ay)  e  ^   - 
n  ^  rv^      TT        n   _,   _^ 


And  then  we  will  designate  the  Dirac  density 
corresponding  to  the  complete  set  of  states  in  the  first 
N  "shells"  as 


N 


N-1 


(II. 1.8)      p''(x,y)  =ZZ  pJ^^y)  = 


n=0 


n 


2  D/2  NJ,  j^  «  ([xl2  +  Iyl2) 

(^)    XZ  ir(ax,ay)  e  ^   - 


IT 


n=0 


n 


where  each  "shell"  corresponds  to  the  complete  set  of 
states  corresponding  to  one  energy  eigenvalue. 
Now 


(II. 1.9)  ZZ  H^(x,y)z''  =  d-z'')    exp 
n=0  " 


o     r  2^  2x  2 
2xy2-(x  +y  )z 


>-T   i/ 


1-z 


Thus 


(II. 1.10)  y 


OO       CD 

z 


m-|_=0  m2=0 


00   m-i  +mo+.  .  Hi-T^ 

ZZ  2 

^D=0 


Hm^K'^l^  «m2^^2'y2^  •••  "mj^K'^D) 


^  Bateman  Manuscript  Project,  Higher  Transcendental 
Functions.  Vol.  II  (New  York,  McGraw-Hill  Book  Co.,  Inc. 
1953)>  p.  194,  eq.  (22) . 


3J> 


2   -D/2 
=    (1-z    )  exp 


2.1     i2,    2-1 


"2x.yz   -    ( |x|    +|y|    )z 


1-z 


or 


(II. 1.11)     ^  z"^  H^^x^y)  -  (1-2^) 

,r>.  ill  r\_,  f\j 


-D/2 


exp 


m=0 


2x'yz-(x   +ly|    )z 


~  ;;.        ~ 


1-z. 


Let 


2        2 


x  =  R+r,        x«y  =  R-r  where  R   =    |R  | ; 


f\,  t\,  r\_,  r\j 


J  =  R   -  V   ,       |x|^+|y|^=2R^+2r^     and   r  =    Ir|> 
«.       "^^       ~  '\'  n.  '^ 


then: 


00 


(II. 1.15)     XZ  z"'H^(R+r,R-r) 
m=0 


f]^      rv_,      ro      «X/      --O 


2     -^/^  r     ^ 

(1-Z   )  expL  J 


2  2 

2Rfz    _   2r   2 

+z  1-z 


-1 


Multiplying  both  sides  by    (1-z)        we  get: 


(II. 1.14)     y~  z^--^  y—  H^(R+r,R-r) 


k=l  m=0 


m  •>-•  '\'  '\'  "v/ 


-1  2   -^/2 

(1-z^      (1-z    )  exp 


2  2 

2Rfz    _   2r  2 

1+z  1-z 


Combining    (II. 1.7)    with    (II. 1.13)    we  get 


OD 


2  D/2  ^   -D/2 


(II. 1.15)      7Z  z'^Pr.fR+^^R-^)    =    (— )         (1-2^) 
m=0 


exp 


1+z       ~       1-2 
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where  p  (x,y)  Is  the  Dlrac  density  of  the  (m+l)-th  shell. 
Combining  (II. 1.8)  with  (II.1.14) 


0°   M  1  M  2  D/2      ,     ^  -D/2 

(II. 1.16)   21:  z    P  (R+r,R-r)  =^  (^)    (l-z)-^(l-z2) 


•exp 


-g^R^d-z)   g^d+z) 


1  +  2 


1-Z 


where  p  (x,y)  Is  the  Dlrac  density  for  the  first  N  "shells." 
Prom  (A.  1.12) 


^D 


(II. 1.17)     P(P,R)   =^      dr  p(R+r,R-r)  e 


h 


21r.p 


O.^       '^^   '"O   'X.I   'Vy 


It   can  be   shown   that 
^D 


(II. 1.18) 


D 


^^^  £  2  2,1+z 

dr   e         ^        e  "^   ^ 


p     r^-zi 

,D     D/2      T    ^  D/2      "  ^2   2  4  +  z^ 


,  D        D 
h        g 


From   (II. 1.15),    (II. 1.1?)    and    (II.1.18)    we  get 


00 


.D 


-D 


(II. 1.19)     ZZ   z^'p^Cp.?.)    =%   (1+z)      exp 


m=0 


m  A,  ~ 


h 


/    2^2^     P_\  rl-2^ 
-(g  R   +-f-2)(lT5) 
fi  g 


Similarly   from   (II.I.16),    (II.I.17)    and    (II.I.18)    we  get: 


00 


(II. 1.20)      ZI   z^-V(P.H) 
k=l 


^   (l-z)-l(l+z)-^. 


exp 


2 

) 

-n  g 


/    2p2,   _9l_\  /l-z^ 
-(g  R   +  ^g-j^^Ti^) 
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4  /    2 

Now  a     =  mk/-ft   ,    so   that: 


2^2 


(II. 1.21)   -g-2   +  a'^R 
■fi  a 


^Ico        2m        2 


c       ~ 


where  co     =  Vk/m. 
c 

Thus    from    (II.I.I9)    and    (II. 1.21)    we   have 


00 


(II. 1.22)     2Z  2"V(p,R) 
iFO        "^  -"  ^ 


7D   ^i+-^'''^^Pf:i^H^P'^^(7TT)l 
h  c 


,    -2H    V 
exp    (  ^jp   )    , 


but 


-D 
(II. 1.23)       (1+z)      exp 


^  H(p,R)(-^) 


OD 


=     ZZ    (--)""  L^"'    (:^H(p,r)), 
m=0  ^        ^^c 


7/ 


where  L  Is  an  associated   Laguerre  polynomial. 

Thus 


pD  _  -2H(p,R)        D-1      u 

(^^•1-24)      PJP,R)    =  ^   (-1)"  exp    (-^^)    \      (4;  ^f£'5^^> 


where  p  (p,R)  Is  the  Wlgner-Moyal   density    for  the 
m  t\,  f\j 

(m+D-th  "shell." 




—'        Bateman  Project,  Higher  Transcendental  Functions, 
Vol.  II,  ...,  p.  189,  eqT  (17) . 
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From  (II. 1.20)  and  (II. 1.21)  it  is  clear  that 


d 


GO 


N-1   N; 


(II. 1.25)  ^ZZ  z    P"(P,R)  = 


dH 


N=l 


^   ^D      -(D+1) 
c   h 


•exp 


^H(p,R)(^) 


,-2H, 
exp(^^) 

c 


but  from  (II. 1.23) 


(II. 1.26)  j^   p^(p,R) 


-    pD     N      -2H(p,R)    D     k    ,    _ 
(^)  ^  (-1)   exp  (   ^-    )    V,(^  H(^,r)) 
c  h  c  c 


Let  u  =  ^-^  H,   dH  =  -tt^  du .   Then 


(II. 1.27)     P^(P,R)  =  -^-4 ^ 


GO 


-u/2 


D 


h 


e  "'   Lj^_-|^(u)  du 


■fnir 


H 


To  get  p(R)  and  t(R)  we  note  that 


00 


D-1 


(II. 1.28)    p(R)  ^  %    f    P  P^^P'S^  ^P 


0 


GO 


N, 


^  J  p  — w 

0 


dp  ; 


and  similarly 
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00 


(II. 1.29)       t(R)    -O^     fp^~^^2^^    P^^l'V    ^P 

0 

ri,  ??      D+2   ^p^(p,R) 

=      -    2m (D+2)     J      P  55 ^P 

0 

but   from    (IT. 1.2?) 

:i  N,      „x  T^  -  ^-  H(p,R) 

(II. 1.50)     -^^ =    (-1)     -p     i^jn^  e 

^  he 


The   corresponding  approximation  In  Part   I  was 

Sp(p,R)           p               p^                          2m(H(p,R)-E) 
(1. 5. 26)  .;  '^     =  ^     ^P     p/^       A.  ( '^^9/-,   )• 


To  compare  (II.I.30)  and  (I. 5. 26),  we  proceed 
as  follows: 

For  N-th  shell,  E^^  =  (N-1+  D/2)ttja^,  so  let 

N 
(II. 1.31)    E  =  (N  -  1  +1  +-|)toc   =  ^to^ 


(II.  1.32)     2mV  =   2m(V-E)  +-7^  (mfiCD  ). 

2       c 


But   since  V  =7^  klxl        and  cd     =  Vk/m   , 
c:         ~  C 


Ivvim               (mh    |VV|) 
(II. 1.35)  2mV=-^^^i^ =     ^ 


(mft  o)p)' 
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Equating    (II.I.52)    and    (II. 1.33)    we  get 


N, 


D 


(II. 1.34)      2^    (mto^)      -    (2m(E-V))(mftjo    )    -[mft|vv|]      =  0. 


This  may  be  regarded  as   a   cubic   equation   In    (mttjo   ) 
The  real   solution  to    (II. 1.34)    Is 

1/3   r         I       i-,2/3 
2.^' -^   [mh    |VV|] 


(II. 1.35)  mto 


N 


iTJ" 


f(a)    , 


D 


where 


(II. 1.36)  f(a) 


1 

;t75 


1   -,_8a£  ^     A   +  1^ 


27N. 


D 


27N. 


D 


+    ^/l    + 


8a' 


1   + 


16a- 


27N. 


D 


27N. 


D 


and 


(II. 1.37) 


a      = 


2m(E-  V(r)) 

wv 

(mh    |VV|  ) 


Now   H  =   H-E   +  N^tioD  /\ ,    so   that 

D      c 


(II. 1.38) 


4h 


fiu 


N      +MHzEl 

D  mfio) 

c 


and   from   (II. 1.35)    we  get; 
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.  2m(H-E)  1/3 

(II. 1.39)    -4^  =   N   +  ■ -27^ (4  N  ) 

^c      ^    (mftlvvD^/^  f(a)      ^ 

but  for  large  values  of  N^^  =  4(N-l)  +  2D  +  2  and  small 
values  of  s: 

8/ 


-u/2  D 
(II.l.HO)   e     Lfj-l^"^ 


and 


(-1)^-^2-^2^^ 


ATF 


Nj^+(%)l/5s         K'' 


A,(s), 


(II.1.41)  f(a)  ^  1   for  large  N^^ 


Thus  from  (II.I.30)  and  (II.I.35) 

rxT  1  U?^     ^P  ^£-S^      -1       2p      „  ,    2m(H-E)   . 
(^^•l-^^)     5^  -  ^   [^IVvll^/^'i^  (m.|VV|)2/3' 


which  aside  from  a  factor  of  2  is  exactly  the  approximation 
(1.5.26).   The  factor  of  2  is  caused  by  the  fact  that  we 
included  a  spin  factor  of  2  in  all  equations  of  Part  I 
but  not  Part  II. 


-^        Bateman  Project,  Higher  Transcendental  Functions,  II, 
p.  200,  equations  (10-(13).   It  should  , be  noted  that  the 
Airy  function  in  Bateman,  A(t)  =  7r/3^/5  A.(-t/3-^/^)  where 
Aj_(t)  is  the  Airy  function  of  the  Handbook  of  Math.  Functions 
used  in  this  paper. 
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To  justify  the  substitution  of  approximation 
(11.1.^2)  in  Integrals  (II. 1.28)  and  (II. 1.29)  to 
compute  p(r)  and  t(r)   respectively,  we  should  note 
that 


P 


D-1 


.z2£ 


A,(- 


2m(H-E) 


2^[mfi|VV|]2/5        1    [mfi|VV|]2/5 


■) 


as  an  Integrand  makes  the  larges  contribution  to 
Integrals  (II. 1.28)  and  (II. 1.29)  in  the  part  of  Its 

domain  where  It  Is  most  accurate  as  an  approximation 

D-1  ^P^^'V 
of  P    ;v- — -   namely  for  H(p,R)  near  E. 

op  ^v  fXj 

To  summarize,  the  exact  Integrals  are 

2m(E-H) 
(mil|VVl)2/5f.(c,) 

(II. 1.45)     pN(p,R)=-%      '^  r       (_i)N+12D^!d//^ 


-CD 


-u/2   D 


Vl^^^ 


u=Nj5-(lfNj^)^/^t 


dt 


(II. 1.44)   p(R)  = 


Q 


D 


h^D 


[m^^[VV| 


D/3 


g 

TTST 


/  tf^ 


D/2      D/2 
t]    [f(a)] 


-co 


(-1)    2  (— )    e    L,,_,(u) 


3 


^N-r-'      dt 
u=N^-(4Np)l/5t 
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(II. 1.45)   t(R)  = 


D+2 
Oj^[mh|VVl]  ^ 

2m  h^(D+2) 


a 


TJ^ 


I    ^f^ 


D+2 


D+2 


t]  2   [f(a)]  2 


-oo 


N+1  D  N^  1/3  -u/2  D 
•    (-1)    2  C-^)    e    Ln-1^^^ 


dt 


u=Nj^-(4Nj^)^/^t 


where 


f(a) 


1 


1  +  -^  +  /I  + 


i6c^ 


3a' 


27N 


D 


27N. 


D 


.  ^/. .  ^' 


27N. 


1  + 


l6a^ 


D 


27N. 


D 


and  as  In  Part  I  a  = 


2mfE-V(r)) 


(n*  |VV|) 


27T 


but  for  large  values  of  N^^  =  4N  +  2D  -2, 


f(a)  ~  1 


N+1  D  N^  1/3   -u/2  D 
(-1)    2  (^)     e    L^_i(u) 


-2/3 
-v  A^(-t)  +  0(Nj^    ) 


u=Nj^-(4Nj^)^/5t 
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Section  2.   Dlrac  Density  for  Harmonic  Oscillator 
with  D  =  2k+l. 

When  D  Is  odd.  It  Is  possible  to  get  an  explicit 
closed  form  for  the  Dlrac  density. 

To  do  this  we  refer  back  to  Part  II,  Section  1. 

2  D/2  D  -a^(R^+r2) 

(11. 1.7)  p^(R+r,R-r)  =  {--)  H^(a  (R+£)  ,a  (g-r )  )e 

(density  of  (n+l)-th  "shell") 

M  2  D/2  N::!   ^  -a^(R^ 

(11. 1.8)  p^(R+r,R-r)  =  (^)    >   H;(a  (R+r )  ,a  (R-r)  )e 


n=0 


and 


oo   £-1   £-1     D 
(II. 1.14)  y         z    y~~   H  (a(R+r),a(R-.r)) 


=  (1-z)   (1-z  )     exp  [ J 

1+z      1-z 


For  D  odd,  D  =  2k+l,  we  can  write: 


00     D 
(II. 2.1)  XZ   z  H  (a(R+r),a(R-r)) 
m=0 


(-1)^   ^   d   vk,   d   ^k,-■  ^2,  ^^^...,.r2a^R^a^   2a^r^z. 
(2a^z)2^  d(R^)    d(r^)  ^"^^     ^"^ 
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(2a^z)^^  m=0     d(R^)    d(r^) 


H^(a(R+r),a(R-r))  . 


Thus 


(II. 2. 2)  H^^+l(a(R+r),cx(R-r))  = 


/  -,  vk    -,   k    ,   k  1 

^-^Mtt  (— V^  ^-^)   "     (a(R+r),a(R-r)) 

(2a^)^^  d(R^)    d(r^)    m+2k 


Also 

00   l-\   &-1      D 


2_  z    > H^(a(R+r),a(R-r)) 


I^      m=0 


2  >2k+l^  w^2v  ^  ^  w  2  J    ^-^'2  J 


(2a ^z)^^^-"  d(R^)    d(r^) 

2   2  2   2 

rTTo^"\                             r2aRz        2arZn 
(II. 2. 3)  exp    [   --^^ 3^3^   ] 


2    ,2k+l''  wr.2^  ■'     ^  w    2-  ■' 


(2a   z)^"^"^^   d(R    )        d(r^) 


CO  I 

y~  z""  H^(a(R+r],a(R-r)). 

m=0  "^ 


Thus: 

£-1     D 
(II. 2. 4)     2Z  H^(a(R+r),a(R-r)) 
m=0 


(2a2)2^+l      d(R2)^      d(r2) 


^^^^   ^:;7^^    ^77^27)  ^l^2^.^^i^^-),aiR-r)) 


kk 


Thus  Pn(R+r;,R-r),  the  Dirac  density  for  the  (n+l)-th 
"shell"  (that  is  the  sum  over  the  ('^p^T^)  degenerate 
energy  states)  for  D  odd  is 


D-1 


(II. 2. 5)   p  (R+r,R-r)  =  {^)  -i^lli e"*^  ^^  +^^ 

(2a  ) 

.  i^L_^(_^__  rr  Vd-i^°^^R+^)^h^+d-i^"^R-^)^ 

d(R2)      d(r2)  2'^+^-l(n+D-l)I 


and 

D-1 

(II. 2. 6)   pN(R+r,R-r)  =  (°£)''^^  Izili.  e'^^^R^^^^) 

(2a  ) 

.  i^^rr,  d   ^^  %+d-i^°^^^+^))Vd-i^"^^-^^^ 

d(R^)     d(r2)  2N+^-^(N+D-l)I 


w 


here  a  =  ymk/h 


N  N-1 

and   p  (R+r,R-r)  =  y~  p^(R+r,R-r) 

n=0 
is  the  total  density  of  the  first  N  "shells"  or,  if  you 

prefer,  the  density  of  the  (^^^"^)  lowest  states. 


^5 


Section  J).      The  Interacting  Harmonic  Oscillator. 

In  this  section  we  will  show  how  a  closed  form  ■ 
for  a  one  particle  density  for  a  system  of  Interacting 
harmonic  oscillators  may  be  obtained  by  a  transformation 
from  the  one  particle  density  of  a  corresponding  system 
of  non-lnteractlng  oscillators. 

The  Hamltonian  for  the  Interacting  system  which 
we  will  discuss  Is 


(II. 3.1)  H 


-  I  p.  I     (mcD  -N  k)  _£  ,   ,p    ^^ ,2 


pj£^    UiicD^^^   P     ,2^k^ 


1=1   ^'"         ^      1=1    -^      '  l,j 
or  equlvalently 


X.  -X  . 


1=1  1=1         1,J 


The  "corresponding"  non-lnteractlng  system  has 
Hamlltonlan: 


N   I   |2     2   ^ 


where  p.  and  x.  are  vectors  of  dimension  D  .    N  Is 
aCI     ~i  P 

the  number  of  particles  In  the  system. 

Transforming  to  center  of  mass  coordinates: 
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N  N 


~    V^r     1=1  -^^  yir     1=1  ~^ 

-        p  _        X 


p.    =      p.    + ^^—  X.    =      X.    + 

^  p  p 


then    (II. 3.1)    becomes 

(II. 3.1.)  H=  (  ^  i|iL  +  i!^  2:f  ixj2 ) 


N.     I-    |2  2      N 


1=1      2m  2     ^    '^1 

2  2 


|pr        mco^-N  k 
^   ^      2m     ^  2  IS'       ^    ' 


and    (II. 3-3)    becomes 

N      I  -    [2  2     N 

(11.3.5)  Ho  =  <  2=  4|-  +  T-  ^  ISil'  ' 

,    ,     l£o'       ^  mo)      I       |2    s 


Now  it  becomes  clear  that  H  and  Hq  would  have  the 
same  wave  functions  if 


|p|2   mcD2_j^^         /ma)2-N  k    IpqI^   mm^ 


2m  ^    2    1^'  '  "  /     2    ^    2m   '   2   '^0 


It   can  be   shown   that 

■Wnk 


> 


'\j 


mcD^ 


^7 


where  U  Is  a  unitary  transformation  written  In  "normally" 
ordered  form: 


D/2     _.  u    /     WIq 

(II. 3. 7)   U  =  p    exp  ^  [1-Plx-  P  .   P  =  7/1-  -^ 


2 
mu) 


Note 


and 


D/2 
U  f(x)   =  p    f(p  x) 


(II. 3. 8)   (f,g)  =  /f(x)  g*(x)  dx  =p^Jf(px)  g*(px)dx 

=   (Uf,Ug)  . 

Written  In  alternate  form  (II. 3. 8)  becomes 


1/2  D  r   t4^  s-t  is-x  It-p 


where  p  =  y  1-  ^^ 

Now  since  any  component  of  p  or  x  commutes  with  any 
component  of  p.  or  x . ,  It  becomes  clear  from  (II. 3-6)  that 


N_  ir  |2     2  N, 


(11. 5.10)  UHpU-l.  (Xf  ^H-JSf-Zf  |x,|2  ) 


1=1  '"'      '   1=1 

T     I  p  I     moD  -N  k    ,  p 

^ .2  ^   2m   +    2     '^'   ^ 
P 
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By  comparison  of  (II.3.IO)  and  (11.3-^),  we  see  that 
although  H  and  UHqU"^  are  not  the  same,  they  must  have  the 
same  elgenfunctlons . 


But  Hq  defined  by  (II. 3-5)  has  elgenfunctlons  of  the 


form: 


(II. 3. 11)  ^i^i'^2' •  •  ' '^d''    ^ 

a       1=1  J=l   [2  ^'J(k.  .)!] 


where  K    (ax..)  are  the  Hermit  e  polynomials  used 

1 J  J    """  "^ 
In  the  last  section,  x  .  Is  the  j-th  component  of  x^ 

and  a^  =  -^  •   For  ease  of  writing  we  will  define  a 

new  set  of  polynomials  in  D  variables. 


D     «n/«  ^P 


(II. 3. 12)  h^(y)   =  TT  p^ 


172- 


It  is  here  understood  that  these  polynomials  are 
ordered  in  some  semi-arbitrary  fashion  ( h^(y) ,h2(y) ,h^(y) . . .} 

such  that: 

First:    The  indicated  sequence  exhausts  all  polynomials 

of  the  form  (II. 3 -12)  and 

Second:   Referring  to  the  notation  of  (II. 3-12) 

D  D 

y~  (n+i).  >  zz  n  . 
j=i       ^      j=i    ^ 


^9 


Now  armed  with  this  new  notation  we  can  write  down 
an  antisymmetric  wave  function  for  E^   of  lowest  energy. 


(II. 3. 13)  ^o(x^,...,x^  ) 

P 


yw 


V      1 


N 


a 


z 


1=1 


X. 


h^(x^),  hgCx^),...,  hj^(x^) 


'^ife) 


^i^^N^' 


\^^N^ 


But  since  H  has  the  same  functions  as  UH„U~  and 
since  the  operator  U  does  not  disturb  the  antisymmetry 
of  the  wave  function,  then  the  wave  function  ^  for  H 
must  be: 


(11.3.1^) 


f     =     JJ   t 


0 


The  corresponding  one  particle  Dlrac  density  Is 


(II. 3. 15)   P(xi,x^)  =  Np  j  ^(Xi,X2,x^,...,Xj^)^(x3_,X2,..x^) 


dXo  dx-.  . . .  dx,, . 
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Combining    (II. 5-9),    (II. 3.13)    and    (II.3.I5)    we  get 

^  ,1/2   2D     N  ^ 

(II. 3. 16)   p(x,,x,)   =   i-^^)       r^(f2)' 

'"^  pa 


rdsdtds'dt'<  Yi  f  yz   f        e"^^ 

ll-*^N   ^1*-^N   ^1"%   ^I'-^N 
Is'x,   It.p     2    ^ 

^  k=l  k 

r 

is  '  'X,   it '  -p  '     o 
k'=l  Jk'  "-^ 


•  5(xo-x^)  5(x,-x^)  ...  5(x,,-Xi,J  dXo  dx^  ...  dx,.  dx,. 


but  '  2.    ^  ■  : 

it  .p 

e^     f(x)  =  f(x  +^^) 

so 

1/2  2D    ^ 
(II. 3.17)   P(x,x')  =  (^^^)   (^)2 

p               74^(s-t-s'.f) 
•/  ds  dt  ds'  df   e^^-P^ "-   I(x,x':s,t,s',t') 


^^  (s-t-s'-f) 

/V   'X/   "Aj     'X/ 
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where 


(II. 3. 18)  I(x^,x^:s,t,s',t')  =^j^jTJ  [dxgdxgdx^dx^.-dx^d: 


5(X2-X2)5(x^-X^)  ...  5(xj^-Xj^) 

Is'x,      2     -fit  o 

(e  '     e        ^  h.  (x,  +  — j  ) 

^k  ^^     /N 


k'=l 


-Is'-x'      2   ,   -nt'  ^ 

Jk'  "^    /N 


-fit 


Now  it  is  clear  that  by  adding  and  subtracting  columns 


1  x^  x^ 
^        2 

1   Xp   Xp 


N 


X-, 


1   X 


•N  ^ 


N 


N 


1   (x-j^-a)   (x-j^-a)' 


1 

• 

(Xg- 

-a) 

(x^- 

a) 

• 
• 

1 

(v 

■a) 

c^- 

a) 

(x-^-a) 
(xg-a) 


.   (xj^-a) 


N 


N 


N 


It  is  also  clear  that  this  can  be  generalized  to 
Hermite  polynomials,  and  then  generalized  again  to  the 


h  (x)  polynomials.   Thus: 
n  '>-' 

Is'x  2  'nt 

N  ^  2    l^k 

A  e/^        e 


tit 


(11.5.19)      ZZ       £ 


^  ^.  ^i  n 


i^  . . .   i,,       i     k=l 
1  N     ^  2  fit 

N       17n~     "    2-'-k+ 
e 


^k  "^     /N 


il---iN       ^ 


^  Vs><%-^'*^^'- ^=-^' 


Continuing  along  this  line  and  using  the  fact  that 

2  -tit   o      is-x,  2  f^t'         "^■^'  '^k 

( II .  3 .  20 )       ( e  ^  e  "^         ^  •  ^  -  '^ 


)-(e 


/N 


o  ^  IS 


i^  '^  I  2 


k'    2  /n''^''^        "^    '    "^ 


'      a  -n 


T^f-|t+t  '-   — (s-s'  ) 


-aM(|t|2+|t'|2' 


?N 


'  /\y  ■        '  f\t 


•    e 


and  the  fact  that: 


o  JK  is        Is '      o 

.c,2|^+_^_(t+t'-  ^  +^^)|2 


(11.3.21)    /    du   e 


2  Vn 


h.(u+-^ 
1  "-      2  /N 


IS  IS 


in 


t+t  '    -  —  +  )    h.(u+ 


IS  IS  ' 

^^"^^   "  "Ed  ^  rmT^ 


6--.(^) 


2   D/2 


We   get    from    (II.5.I8) 


2^2 


^mzll     a^(N-l)|t+t'-  ^(s-s-)'^ 
(II. 3. 22)    I   =    (^) 


,q^  N        2  ,  'W 


'^'    '^  mCO    '^'    "^ 


-#(-.(i*i^^iri% -"^i^^^;;!'^^;^'^^^^/^'  "^ 


(e 


N 


-tl 


>        h.(x,+  -i^-(t+t 

4 1      ~1  -  /TT     ~      ~ 


t        _ 


1=1 


2  /N 


^  Is      is' 

-^  +^^)  h,(x'  +-^(t+t.'-  :£  +  t;^)) 

TiYn         mo)  1  ~1 


IS  IS 

moo       im) 


2  /N 


'K.  ~       ino     ma> 
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Now  let  us  substitute 


r:/       = 


t  +t' 


</ 


s  -  s 


yt=   t  -  t' 

r\j  t\j  t\j 


Xi-  =   s  +  s 


Combining  (II. 3-22)  and  (II. 3. 1?)  with  the  new 
variables  we  get 


^  ^1/2  2D    D/2  p 
(II.5.23)   P(x,x')   =  (gflflpT^   ^^2^    j^^^^^^^ 


2^2 


2.2,1^1 


«  (£^,.  siV)  ^,N-1)17^4^J=^  -«t^ 


+ 


1:71^) 


TT^ 


•  e  e     P  e  P 


>/N 


Ixi' 


•  (x+x' )  - — ^  '(x-x' ) 
~  ~     2  yN     '^ 


N 


^ 


>   h. (x  + 

fe  ^  -^    VN 


(5^-  -^4))  h.(x-  +-^  (^-  -i-J) 


where 


r-i  4\^    =  (r-;i;</)-(^-^/) 


~   intD 


~   mo)  ~      ~   mo)  ~ 


5^ 


The  variables  -<*'  and  Jt  may  be  Integrated  Immediately 
giving  us : 


(II. 3. 24)   p(x,x') 


D 


r  "^  Q  \  r 

^  1/2  H  ol^   ^""^  ^ 


a 


2 


(1-p^) Ix-x- 


4N 


"d^d^[exp  (-7^^-  -^+^  (N  -1)|/-  -i->^l2 


a 


-h^l^l^  -  ^    (|x|2  +  Ix'l^)  +9^  x7^-(x+x') 

^    P 


p 


N 


p 


zi  h,  (x+  -5z  (^  -  ;^>^) )h,  (x'+  ^^-  A^^) ) 


1=1 


1  '^ 


^  p 


~   moo  ~   1  ~ 


~   mo)  ~ 


Now  substituting: 


X 


R  +  r 


X 


R  -  r 


■fi 


£  (v  -  R) 


0_,       r\_, 


a     vN  u 
p  ~ 


and  using  the  fact  that 

D/2   -|!(|x|^+|xM2)  J!e 


a 


>"  h. (x)  h. (x-) 

■'=: r-    1  A^     1  ~ 

1=1 


we   finally  arrive  at   the   equation; 
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a^(l-p^)|r|2 


(II. 3. 25)  p(R+r,R-r)  = 


r\_,       ^\j       r\_,       r\j 


f— ^^^ r)       e 


N 


P 


2,  ,2  ,  2iN  a  p 


rr  2i  i2     ^i^va  p 

//  du  dv  exp  [ -N  a   u    +   /-.^^  ^   u 


V 


2iN  a^P 


p„(v-iu+r,  v-lu-r) 


'0 


*\j  '\j     ^\j  'V*     o^   o_* 


In  the  special  cases  that  N  =  (■'^~^"^^),  N  =  1,2, 

we  have  a  closed  shell  structure  and  we  may  obtain  a 

more  explicit  form  for  p(R+r,R-r). 

Following  the  notation  of  Section  2,  we  will 

N, 


denote      such  a   density  by  p    (R+r,R-r) 
Now   from   equation    (II.I.16) 


00 


t; —     N-1   N,„  ,      o      ^ 
> z        pp,(R+r,R-r) 

N=l 

2  D/2 
=    (^)         (1-z) 


(D+2) 


-D/2 
(1+z)  exp   [ 


1+z 


2  2,,,     X 
g  r  (l+z)i 

1-z       J 


Thus 


00 


(II. 3. 26)      y~  z^-^   p^(R+r,R-r) 


N=l 


i\,     n_»     ""x^     n^ 


-D/2      Nap      D 
(1+-)  (^f^)      exp    [ 


g^(l-p^)|r 


N 


P 


Po  2        2        ^^^  "   P  2iN  a   p 

JJ   d^  dv  exp   [-N^a    |u|      +  -jj^  u-v-      (^P^)      u  •  R] 


exp 


^  4feS^    Hvl^   -    2iu.v-    |ul2)    - 


a^(l+z) 
(1-z) 


Ir.l'] 
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Carrying   out    the   indicated    Integrations   we  arrive  at 
the   equation: 


00 


(II. 3. 27)     ZZ  z^"^p^(R+r,R-r) 


2  D/2        N     S  D/2 


N=l 


r\j     r\^      f\^     r\f 


(^)  ( 


2      2 


) 


exp   [- 


a^(B^-l)    ,,2 
N 
P 


r'^^]      exp   [-  a^( g-5^-)R^   -  aV    ] 


(1-z) 


'    exp 


2n2  a^p%2        N  fi^^^^.i 
— P^ (    P  2) 


N   R    ^    -1 
1    +    (    P    p      p      )z 
NpP^-p^+1 


o   2      2 
ag    zr 

(1-z) 


-^ 


I 


where  we   have   substituted  R  =    |r |    and  r  =    |r|. 


Now   if  we   substitute  7 


& 


2-1 


NpP 


^  and  expand  the 


part  of  the  last  equation  within  \      \    1 
Laguerre  polynomials  and  then  we  get: 


n   terms   of 


00 


(II. 5. 28)   YZ  z^"^  p^(R+£.R-r)  = 

■  '  f\j  Oj  f\.l  f\.i 


N=l 


f\j      r\j      r\j      r\j 


2  D/2  ^      D/2  .         a!2^      2    ,  .        ^.V  2   2, 

D-2 


[ 


00 


i=0 


D/2 


(2a2r2))    (  ^    (-z)"\l±I)"^  ^T^^^) 
m=0  ^    ^  1-7 
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Thus  for  the  Hamlltonlan, 


H  = 


1=1 


2m 


N 
0  ^-£ 

2 


1=1 


1,J* 


we  have  arrived  at  the  Dlrac  density  p  (R+r,R-r)  where 


2  D/2      -D/2 


N, 


(II. 3. 39)  p"(R+r,R-r)  =  (^)    (I-7) 

•     "-v  ro  n_/  n_,  7/ 


exp  [  -  ^^^^  -  a^(l-(N-l)7)  ^2  , 
exp  L    (1-7)       (1-N7)         ■' 


^-1     .m,l+7- 


D-2 


ZI  (-1)"  (f^) 

m=0        ^   ^ 


7    m    ^ 


-2   ,2a2R^  y^      (2aV) 


1-7 


2^  ^N-l-m^ 


-/■- 


1  -y  1+  ai>^ 


where  7  = 


N 


ko      .2 
— ^=^ —  and  a 


V^ 


0 


-fT 


/ 


-^ 


To  conclude  this  discussion  cf  the  exact  solution 
for  the  Interacting  harmonic  oscillator,  we  will 
consider  the  lowest  energy  level. 

Rewriting  equations  (II. 3. 4)  and  (II. 3. 5) 


N    I-  |2 

-   Ip  ' 


1 '     mo) 


(11.3.4)  H  =  (  Zf   2^ 

1=1 

N   I-  |2      ^   N  I 

r  (  Z_  -^5^  +  ^-  zZ  lxJ  )+(-p- 


(II. 3. 5)  H, 


0 


1=1 


2   1^  -^ 


2m   "^"2   1^0'  ^• 
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Examining  these  two  equations  we  see  that  difference 
in  energies  arises  from  the  second  terms  thus 


E  -  Eq   =  1^  [  /co2  _  ^  -  0.  1 
or 


For  the  closed  shell  system  N  =  (  -^  )  •  It  can 
be  shown  that  after  summing  over  all  degeneracies  that 
for  the  total  energy  E^  for  the  non-interacting  system 
(II. 5. 5)  we  may  write 


D(D-1+2N)     D-l+N 

E„   =   (         )  ^^  CD 

°       2(D+1)       D 


Thus  for  the  Hamiltonlan: 


N    I   i2    ,     N 


^0     ^  |2  ^  k  ^ 


H  =  >     -^ —   +  -^   >    Ix^  r  +  4f  >    |x.-x. 


the   lowest   total    energy 


D(D-1+2N)    ,D-1-^N.^     /VV    ,    D^,    r    PQ        /5^] 
-        2(p+i)        ^D       ^^^        m         ^2"Wm       V         m      ^ 
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Section  h.      The  Hartree-Fock  Solution  Minus  the  Exchange 

Term  for  the  D  Dimensional,  N  Particle 
p 

Interacting  Harmonic  Oscillator. 


Consider  the  Hamlltonlan 


(II. 4.1)   H 


N 


1=1 


2   2 

±L-  S/         +  — 
2m  X.    2 


1     N 
^0  ^^ 


^-  1^1 


+ 


ir  ^^    Ix.  -X  .  I 


i,J 


where  x^  =  ^^n'  ^±2'       "'   ^n) 
,2 


and  V 


2 


^       a2 


~i   m=l   ox. 
im 


2   . 


Let  us  rewrite  (II. 4.1)  In  the  form: 


N 


J2.     2    k^+(N  -l)k  _%     o    

1=1      ~.  1=1        i,j 


Suppose  we  make  the  Hartree-Pock  approximation  that 


^ter^2'"-^p^ 


"^   P 


^2^1^ 


P 


•    ^l^^N    ^ 


P      P 


then 


(II.4.3)<H>=   -liZZ    /<(x)   v2  ^.(x)    dx 

1=1     '^  ~ 

y~  [ i^U^)  ixi^  ^. (x)  dx  - 


ko+(N^)k 
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k>    /  ^*(x)  ^"^CxMx-x'  (^.  (x)^.(x')-^.  (x')^.(x))  dx  dx' 

1  ^  J 


i>j 


Now; 


5«H>+2i:  A^j  j  ^*(x)  ^j(x)  dx  )  =  0  . 

And  since  the  coefficient  of  6^. (x)  is  equal  to  zero. 


we  may  write 


^2   o  l^n+(N^-l)k  ,    o 


-  k  Jp(x-,x')x-x'  dx'  ^^(x)  +  kjrp(x,x-)x-x'  ^.(x')dx' 


Since  (?\.  .)    is  a  Hermitian  matrix,  it  may  be  diagonallzed  by 
a  unitary  transformation  ((J  j^  J   that  will  give  us  a  new 
set  of  orbital  functions 


^i  =  ZI^,kVo-  ^i,J  • 


We  then  get  the  new  set  of  equations 


6l 


^2   o         k^+(N  -l)k    o 


-  k  /  p(x',x')x-x'  dx'  u.(x)  +k  /  p(x,x')x-x'  u.(x')dx' 


=  e.u. (x) 

1  X  ^ 


where  p(x,x')  =  ^^  ^.  (x)  V'^-  (x')  =  >   u  (x)  u.  (x' )  . 

1  1 

However    /  p(x',x')x'X'  dx'  =  0  , 

and  the  last  term  on  the  left  hand  side  of  (II. 4. 5)  Is 
the  contribution  from  the  exchange  term  which  we  here 
choose  to  Ignore. 

We  are  thus  left  with  the  equation 

J2.     r.  k^+(N  -l)k    o 

-|-v2u.(x)  +  Q   P Ixl^u.Cx)   =  e,  u,(x)  . 

2m  X  1  ~        2      '~'    1  ~       i   1  ~ 


We  thus  get  the  usual  eigenvalues 


^  ykQ+(N  -l)k 

^i  =   ^^i  ■"  2^  ^  y — ^ — 

n.+D-l 
with  a  degeneracy  of  (      ) . 

D-1 
If  we  continue  to  ignore  the  exchange  term 

N 

<H>  ^   7^  e.  . 

1=1   ^ 

For  N  "shells"  we  have 
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N-1  n.+D-l  N+D-1 

>  (        ^  )      =      (  )      different   orbltals 

n.=0  D-1  ^ 

1 

thus    If  N     =    (  ~      )    ,      N  =   1,2,...,      we  may  write 

N-1        n.+D-l  ^  /k^+(N   -l)k 

<H>-  II  (  i-i  )  K+f)-^  y— -rf — 

n.=0  ^ 

1 


or 


(11.4.6)    <H>^    "tillT)^     (%'^^)-fi/ 


V^p-i"^ 


m 


For   comparison   let  us   write   the   exact   energy   from  page   59: 


k^+N  k       T^  /k^  /k^+N  k 


/MN          D(D-1+2N)       ,D-1+Nv    .    /^0""p"    ,    D  *.    r    /--2         /    0      P      , 
<H>  =     ~  2(D+1)         ^      D      ^   ^/ m         +  2  ^   ^  V—  -  -/       m  J 


Section  5.   The  Thomas  Fermi  Approximation  for  the 
Interacting  Harmonic  Oscillator. 

The  Wigner-Moyal  density  without  the  spin  factor 
of  2  for  the  Thomas-Fermi  approximation  is: 

-,  2 

P^^'V      =  — D-  e  (E  -  V(r)  -  ^)    . 


Thus  _  ^ 

D+2 


°°  2  D-1       Op   [2m(E-V)]  ^ 


(II. 5.1)  t(r)  =C^/  ^   p(p,r)  p    dp  =-^  ^^ 


0  ^2"^ 
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On  the  other  hand 

CD  ^/2 

D-1       Op   [2m(E-V)]. 


(II. 5. 2)   p(r)  =  r^  /  p(p,r)  p    dp 


b  h 


We  thus  arrive  at  the  approximation 

D+2 


D+2 


(II. 5. 5)   t(r)=|^  §7^ [p(r)]^   . 

(Op)    (D+2) 

The  Hamlltonian  used  In  the  last  section  was 

_%    Ip.  I^    k„  +  (N  -l)k         2 
(II.5.M   H  =  Zf -ij^  +  ^^-^E EZIx^l-k^Zx, 

1=1  1        1,J 

i>J 
Now  using  approximation  (II. 5. 3)  we  find  that  the 
total  energy  of  the  system 

D+2 


X  . 


2      ^  D 


D+2 


+  ^L_E r|xl%(x)  dx 


r\j  '\j 


-  -o  /  x-x'[p(x,x)p(x' ,x' )-p(x,x' )p(x' ,x)  Idx'dx 

Now   /  x-x'p(x,x)  dx  =  0 
and  we  will  Ignore  the  exchange  term  as  we  did  in 
the  last  section.   This  gives  us 


6h 


2     ^  D 


D+2  p^2 


(^^•5-^^<H>=  1^   ,,  "2/D,.,.,  /tp(x)]^   dx 


(Op)  "/"(D+2) 


^kQ+(Np-m 


/  |x|   p(x)  dx  . 


Using  the  relation  that 

5(  <H>  -   E  f  pix)    dx)   =   0  , 
we  obtain  the  equation 

v,2    T,^^^  2/D   k  +(N  -l)k    2 


or 


O^        k^+(N^-l)k    p  D/2 
^  [2in(E-  -^^£ Ixl^)] 


(II. 5. 7)   p(x)  =  -D^ 

h  D 


We  can  obtain  E  from  the  equation 


00      D-1 
Op  /p(x)  X    dx  =  Np 

0 


(Cu)2[2mf/2'X        D/2  D-1 


where  A  = 


V(%-1^^ 


2 
If  we  substitute 


X  =  ./|  sin  0 


we  get 
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^/2        7r/2 
0 


(^p)  [^"^1       E^    r  D+2  .   .  D-1 


We  now  use  the  relation 
■w/2 


'n\    r^  fVay 


J  sin'^-l   9    cos"-l   9   dO      =     I     -  ^^^^ 
0 

^  r(|) 

Prom  these   equations   we   obtain 


1/D                                     /k^CN   -l)k 
(II. 5. 8)      E  =  -^   o)  (DI   Np)  where     03  -    '  ^ 


m 


but  using  this  result  with  (II. 5. 8),  and  carrying  out 
the  Indicated  integralons  of  (II. 5. 6)  we  obtain  the 
equation 

2+1     /k-+(N  -l)k 

,N+D-1n 
For  the  closed  shell  structure  where  N  =  (  j^      ) 

we  get  the  Thomas  Permi  approximation 


D    r(N4-D-l);.  .  .(N+D-1):,  1/"  /■^0^'V^"' 


(II. 5. 9)  <H>  =  -^prpiyy  [\n-1);  ^  ^  (N-D!  ^   ^ 


m 
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And  for  comparison  the  Hartree-Pock  minus  the 
exchange  term  discussion  in  Section  4  gives: 


<H>  =  TDTryy  ^    (n-i)i    ^  '  ^n+  ^-)  -^  / — 


m 
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Appendix  1 

Part  I:  Wigner-Moyal  (Dirac)  Density 

Notation: 

u       scalar  quantity 

u       vector  quantity 

u       operator  with  one  component 

u       operator  with  several  components 

<P(p,x)>  expectation  value  of  operator  P(p,x) 


^/  ~ 


Consider  a  Moyal  ordered  operator  P(p,x) 
(A. 1.1)      P(p,x)  =  ZZ  Pv  m^P  '^^'^ 


where 


sum  of  all  possible  orders 


(A. 1.2)   {p\x"^}  =  [1/(^^)1  •(      _     "    .  ) 

of  k  p  's  and  m  q  's 

If  we  limit  ourselves  for  the  moment  to  an 
operator  for  which  a  fourler  transform  exists,  we 
may  write: 

(A. 1.3)   P(P,x)  =   /TdT  dO  f(T,0)  exp  (irp  +  lOx)  . 

Using  the  corresponding  scalar  function 
(A. 1.4)      P(p,q)   =  ZZ  Pk,m  P   ^ 
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(A. 1.5)   f(T,0)  =  (l/27r)^  //  dp  dq  P(p,q)  exp  (-Irp-lOq) 


so 


(A. 1.6)  P(p,x)  =  (l/2Tr)^  //  dT  dO  dp  dq  P(p,q) 


•  exp  (iT(p-p)  +  10(x-q)) 

(A. 1.7)    P(p,x)    =  (l/27r)^  //  ^^  ^'^   ^®  ^P  ^^  P{p,q)^*(x)- 

exp  {lT(p-p)+10(x-q) )^(x)  . 

Here  we  use  the  equation: 

iTp+lOx      -^  rp   10X  "5  Tp   -2/ 
(A. 1.8)      e         =  e     e    e"^    ; 

(A. 1.9)   <P(p,x)>  =     ■  ^   r,      /Tdx  dT  d9  dp  dq  P(p,q)- 

{2irf    ^ 

-^     ,      -lTp-10(x-q)  ^ 
.   (e   ""  ^  (x))  e  (e  '^  ^(x)). 

Noting  that  e^'^P/^^(x)  =  e^'^^^^  ^/^^^(x)  =  ^{x+ ^x/2) 
and  then  Integrating  out  0  and  then  q,  we  get: 


(A. 1.10)   <P(p,x)>  =  -—  JJ  6x   dp  P(p,x)  J  6r   /(x-^) 

-iTp 


e   ^(x+  -5-)  . 


-^W.  0.  Kermack  and  W.  H.  McCrea,  Proc  Edinburgh  Math.  Soc. 
Vol.  II,  Series  2,  (1931),  P-  224. 
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From  this  we  see  it  is  natural  to  consider  the  phase 
space  distribution  function: 

(A. 1.11)   p(p,x)  =^  J  dx   /(x-^)  ^(x+^)  e 

This  is  generally  known  as  the  Wigner  or  Wigner-Moyal 
density.   Of  course  this  may  be  easily  generalized  to 
multiple  dimensions,  multiple  states,  and  multiple 
particles . 

(A. 1.12)   p(p,x)  =  J^JJ   dr  p(x  +  2^  ,  X  -  ^)e     '"  --, 

where  D  equals  the  number  of  dimensions  and  the  better 

10/ 
known  Dlrac  density. 

Rewriting  (A. 1.10),  we  have 
(A. 1.13)  <P(p,x)>  =   //dp  dx  P(p,x)  p(p,x). 


To  generalize  this  relation  for  operators  with  no 
fourler  transform  we  can  use  the  generating  operator: 


(A. 1.14)   P(p,x)   =   /J  dT  dQ  e 


i(Tap+9bx)   -T^-0^ 


(The  coefficient  a  d"^  will  be  some  constant  multiple 


„  r  -k  -m-i 
of  [p  ,x  } 


^^L.   D.  Landau  and  E.  M.  Llfshltz,   Quantum  Mechanics, 
(Reading,  Mass.,  Addison-Wesley 's,  195»),  p.  35  ff- 
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Thus  we  see  equation  (A.I.I3)  wHl  apply  to  any 

operator  P(p,x)  of  the  form  P(p,x)  =  XZ  Pv  m'-^  '^  ^ 

-V   -m  ^'^        ' 

If  y       p    <[p  ,x  }>  is  convergent. 

k,m   '"^ 

It  Is  now  clear  that  we  may  expand  this  last 

generalization  to  multiple  dimensions  so: 
if 

P(P,x)   =  ZZ  Pk  m^P  '"^   ^ 

k,m  ^'      ^     -^ 

and 

kTS  ^'    ^  ^ 
is  convergent,  then 
(A. 1.15)     <P{p,x)>  =   /Tdp  dx  P(p,x)  p(p,x) 

Where  p(p,x)  is  the  Wigner-Moyal  density  defined  in 
equation  (A. 1.12). 

Note  I   I  have  used  here  the  same  functional  notation 
for  the  Dirac  and  Wigner-Moyal  densities,  but  this  should 
not  create  confusion  since  the  Wigner-Moyal  density  is 
a  function  of  both  momentum  and  coordinates,  while  the 
Dirac  density  is  a  function  of  coordinates  only. 

Historically,  there  is  interesting  comment  to  be 
made.   P. A.M.  Dirac  appears  to  be  actually  the  first 
man  to  use  this  phase  space  density  referring  to  it  as 
a  "pseudo  classical  density"   [l].  Dirac's  "pseudo 
classical  density"  was  not  properly  normalized  and 
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Dlrac  failed  to  observe  that  expectation  values  for 
measurable  quantities  could  be  computed  directly 
without  converting  back  to  the  Dlrac  density. 

In  1932  E.  Wlgner  published  a  brief  paper 
specifically  on  this  density  without  reference  to  the 
Dlrac  paper  [lO],   He  noted  that  this  density  Is  not  a 
genuine  probability  distribution  (it  is  always  real  but 
It  can  take  on  negative  values).   More  Importantly 
however  he  Indicated  that  for  many  computational  purposes 
this  phase  space  density  can  be  treated  as  if  it  were  a 
genuine  probability  density. 

For  several  years  afterwards  this  phase  space 
density  was  known  as  the  Wlgner  density.   Then  in  19^9 
J.  E.  Moyal  published  an  extensive  paper  on  this  phase 
space  density  [11].   (Most  of  the  arithmetic  in  the 
opening  paragraphs  of  this  appendix  were  lifted  directly 
from  Moyal 's  paper.)   Since  the  publication  of  this  paper, 
this  density  has  generally  been  referred  to  in  the 
literature  as  the  Wigner-Moyal  density. 

One  of  the  things  that  Moyal  observed  was  that  this 
density  satisfies  an  equation  very  close  to  the  classical 
Llouville  theorem,  namely: 

•  H(p,q)  P(p,q,t). 
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This  is  a  pretty  equation  and  it  gives  some 
indication  of  how  close  the  Wigner-Moyal  density  is  to 
the  classical  density.   However  as  with  all  such  dynamical 
equations,  it  is  necessary  to  insert  proper  boundary 
conditions.   For  the  corresponding  classical  equation: 

^F  ,   ^  ,  ^    hU        hF  SH  SF 

all  non-negative  functions  F(p,q,0)  which  are  first  order 
dlfferentiable  in  p  and  q  are  physically  admissible  as 
initial  boundaries.   Also  any  non-negative  F(p,q)=f (H(p,q) )  , 
where  f  is  arbitrary,   is  physically  admissible  as  a 
stationary  solution. 

However  in  quantum  mechanics  this  is  no  longer  the 
case.   For  instance  for  the  harmonic  oscillator  f(H(p,q)) 
where  f  is  arbitrary  is  a  stationary  solution  of  (A.I.16) 
but  may  not  correspond  to  any  Wigner-Moyal  density. 

It  is  not  at  all  clear  what  boundary  conditions  are 
necessary  for  (A.I.16)  and  without  them  it  becomes  almost 
useless . 

Despite  this  drawback  several  physicists  In  recent 
years  have  found  the  Wigner-Moyal  density  useful  for  some 
kinds  of  computations,  just  as  Dirac  did  back  in  1950- 

In  1957,  Kompaneets  and  Pavlovskil  published  a  paper 
showing  that  the  Thomas-Fermi  density  for  the  atom  with 
Welsacker  and  exchange  corrections  could  be  extracted  from 
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the  Hartree-Pock  equations  [j>]  .      In  this  paper  Kompaneets 
and  Pavlovskll  referred  to  the  "pseudo  classical  density" 
of  Dlrac's  1950  paper  and  then  rederlved  the  Moyal  bracket 
equation  (A.I.16).   From  an  esthetic  point  of  view  there 
were  two  shortcomings  In  their  paper.   The  first  shortcoming 
was  that  Kompaneets  and  Pavlovskll  were  apparently  unaware 
of  the  work  of  Wlgner  or  Moyal  and  so  there  Is  no  Indication 
of  the  fact  that  It  Is  unnecessary  to  convert  back  to  the 
Dlrac  density  In  order  to  compute  expectation  values.   The 
second  and  perhaps  more  serious  shortcoming  resulted  from 
their  use  of  the  Moyal  bracket  equation  (A.I.16).   They 
"solved"  the  equation  without  Introducing  the  necessary 
boundary  conditions  and  thus,  they  were  unable  to  show 
that  their  solution  was  necessarily  the  right  solution. 

In  i960,  Baraff  and  Borowltz  corrected  the  second 
shortcoming  by  using  what  amounts  to  a  Wlgner-Moyal 
Green's  function  [9].   However,  Baraff  and  Borowltz  also 
did  not  relate  their  work  to  that  of  Wlgner  and  Moyal 
so  their  paper  also  suffered  slightly  from  the  first 
shortcoming  In  the  Kompaneets  and  Pavlovskll  paper. 

In  closing  let  us  look  at  the  Wlgner-Moyal  densities 
for  a  few  simple  problems. 


For  the  minimum  wave  packet: 

r    (x-<x;     

2  "   -fi 


(A. 1.17)     ^(x)  =  [27r(Ax)2]'    exp  [  -  ^^"^^^j>   +  4§^  ]~ 

4(Ax) 


— ^Schiff,  Quantum  Mech. ,  p.  57 
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Using  equation  (A.l.ll)  we  compute  the  corresponding 
Wlgner-Moyal  density  and  get: 


(A. 1.18)    p(p,x)  =  ^  exp  [  -  ^^-<^>/   -    iP-<P»^   ] 

2(Ax)        2(Ap)^ 


irfi 


where 


and 


(Ax)^  =  <'(x-<x>)^> 


(Ap) 


2 


<(P  -<P>)"> 


,.    ^2,.    ^2         ^ 

( Ax )  ( Ap ) 


T" 


In  the  cases  of  the  free  particle  and  the  harmonic 
oscillator,  the  Moyal  brackets  reduce  to  the  classical 
Polsson  brackets.   Thus  for  these  two  cases  we  should 
expect  some  very  classical  looking  solutions.   Using 
equations  (A.l.ll)  and  (A. 1.12)  and  following  the 
techniques  used  by  Schiff  on  pp.  56  to  58  and  pp.  67 
to  69,  we  get  the  time  development  of  the  minimum  wave 
packet  for  the  free  particle  (one  dimension) 


(A. 1.19)      p(p,x)    =  -^  exp   [    - 


(P-^P0»^        (x-<Xo>-  ^)2 


irfi 


where 


<Po> = <P> 

^  '  t=0 


0 
2(APq)' 


<Xq>   =    <x> 


0" 
2(AXq)' 


] 


t=0 


(APq)2   =    <(p-<p>)2> 


t=0 


,    (Ax^)^   =  <(x-<x»2> 


^0 


t=0 


f\ 


(APq)"(AXo)^   =-^   , 
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and   for   the   harmonic   oscillator    (one   dimension) 

(A. 1.20)  p(p,x)    =  -Jc  exp   [    -  l2^^<2^  -  i£z<£>l!   j 

^^  2(Ax)2  2(Ap)^ 

where 

<Po> 

/  x>  =  <'x^>   cos  cot   H sin  cot 

^  ^0  miD 

k<'xQ> 
<p>  =  ^Pn'^    cos   cut    -  sin  aot 

and  where  cd  is   the   classical    frequency  =  /k/m   . 

In  his  textbook,  Schiff  mentions  the  fact  that  in 
coordinate  space  the  free  particle  minimum  wave  packet 
spreads  as  time  increases . 


(A. 1.21)    p(x)  = 


P   (APo)2t2 
2Tr({AXo)2+  % ) 

m 


<Po>t  2 
exp  [ g  g   1  . 

p       (APf^)   t 

((Ax  )2  +  ^2—) 

m 


However  since  the  Moyal  brackets  in  this  case  collapse 
to  Poisson  brackets,  the  Wigner-Moyal  density  satisfies 
the  conditions  of  the  classical  Liouvllle  theorem.   Thus 
as  is  apparent  from  the  form  of  (A.I.16),  the  minimum  free 
particle  Wlgner-Moyal  wave  packet  moves  along  the  classical 
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path  in  phase  space  and  it  "spreads"  only  in  the  sense 
that  It  changes  shape. 

Referring  to  (A.I.16)  and  (A.I.19)  let  us  consider 
the  ellipse: 


(P-<P0>)2    ^   (x-<x,>-f)2 


2(APq)^  2(AXq)^ 


In  the  figure  on  the  following  page,  since  the 
area  of  the  little  parallelograms  remains  invariant 
as  t  increases,  we  see  by  an  obvious  limiting  process 
that  the  area  of  the  ellipse  remains  invariant  as  t 
increases . 
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Part  II:   Klrkwood  Density  [12]. 

The  arithmetic  for  the  Klrkwood  density  Is  much  the 
same  as  that  for  the  Wlgner-Moyal  density.   The  chief 
difference  Is  that  for  Klrkwood  computations  one  begins 
with  a  "normally"  ordered  operator. 

(A. 1.22)     P(p,x)  =  ZZ  Pv  r.  ^^  P""   • 

k,m  ^'^ 

Assuming   for   the  moment   that   P(p,x)    has   a   Fourier 
transform  we  may  write 

CT  i®^      i''^P 

(A. 1.23)  P(p,x)    =      //   dT  dO   f(T,0)    e  e 


Now  using  the  corresponding  scalar  function 
P(p,<,)  =^P^,^q%" 

we  get: 

rj^  10q   Ixp 

(A. 1.24)     P(p,q)   =   //  dT  dO  f(T,9)  e    e 

-,      rj"  -lTp-10q 

(A. 1.25)      f(T,0)   =  — ^   //  P(p,q)  e         dp  dq 

(27r)^  ^J 

1         rr  i0(x-q) 

(A. 1.26)      P(p,x)   =  — ^   //  dT  dO  dp  dq  e 

(27r)^  JJ 

It(p-p) 
•  e        P(p,q) 
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(A. 1.27)    <P(p,x)>  =  — ^   //  dT  dO  dp  dq  dx 

(27r)^  uLi 

^     19(x-q)   It(p-p) 
•  P(p,q)  ^  (x)  e        e       ^(x) 


(A. 1.28)    <P(p,x)>   =  — ^   //dT  d9  dp  dq  dx 

(27r)^  ^' 

^  10(x-q) 

•  P(p,q)  i/    (x)  e       ^(x+fir) 


Integrating  out  0  and  then  q  we  get 


(A. 1.29)    <P(p,x)>  =  -^   /7  dp  dx  P(p,x) 


■        6x  Tj/   (x)  e     ^(x+hr)  . 
We  thus  define  the  Klrkwood  density  Pj^(p,q) 

+0O 

(A. 1.50)     p^(p,x)  =     ^    J      ^  (x)  e     ^(x+fir)  dx 

-00 

Then  from  (A.I.29)  and  (A.I.30) 


(A. 1.31)    <P(p,x)>  -   /Tdp  dx  P(p,x)  Pj^(p,x) 


In  the  same  fashion  as  we  did  for  the  Wlgner-Moyal 
density,  It  Is  possible  to  extend  equation  (A.I.3I)  to 
apply  to  operators  for  which  no  Fourier  transform  exist 
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And  as  we  did  for  the  Wlgner-Moyal  density,  we  may 
generalize  the  Klrkwood  density  to  multiple  dimensions, 
multiple  particles,  and  multiple  states. 

Last  we  should  note  two  equations  used  in  the  main 
body  of  the  paper  (A.I.5J1)  and  (A.I.38). 

First  taking  the  Fourier  Inverse  of  (A.I.50)  we  get; 

00         .  , 

(A. 1.32)     J   pj^(p,x)  e     dp  =  ^  (x)  ^(x+1Tt') 
-00 

Substituting  x'  =  x+hr '   we  have: 

(A. 1.33)     f   (x)  ^(x')   =   /  p^(p,x)  e  dp 

-00 

Taking  complex  conjugate  of  both  sides 


00  1  p^y'^ 

(A. 1.34)     ^(x)/(x')   =  r    pj^(p,x)  e  dp, 

-OD 

Second,  to  get  the  operator  which  transforms  the 
Klrkwood  density  into  the  Wlgner-Moyal  density 


-i-ti  ^^  -ifi  h^        ^ 

2  "Sp^  ^  2  ^p^   ® 


(A. 1.35)      e         pj^(p,x)  = 


e 


2Tr     J 
-00 


dT 


-irp  ^ 
e     ip   (x)  ^(x-rtfr) 
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-±n    y 


+CD 


(A. 1.56)      e    2   ^^^  pj^(p,x) 


27r 


dr  e 


-00 


^    (x)    ^(x+^t) 


-Ifi   s^ 


00 


-iTp 


(A. 1.37)      e    2    ^^p^(p,x)    =^/dT   e"  /(x-^)^(x+^) 


Thus 


(A. 1.58) 


-iff      a^ 

2    6pdx        /         \ 
e  ^        p^(p,x) 


=      p(p,x) 
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Appendix  2 

Here  we  have  essentially  copied  from  the  paper  by 
W.  Macke  and  P.  Rennert  [6]. 
Consider  the  operator 

-uA    uA 
(A. 2.1)      F(u)   =   e    Be 

where  A  and  B  are  operators  and  u  is  a  scalar. 


P(0)  -  B 

Oil* 


Now 


-uA   uA    -uA   uA 
Ft(u)  =  -  Ae    Be   +  e    Be  A  =  [F(u),A] 


(A. 2. 3)       F"(u)  =   [P'(u),A]   =   [[F(u),A],A] 
F^5^(u)=  [[[F(u),A],A],A] 


.(n) 


[1   2   .  .  .  n  1 


n. 


(A. 2. 4)       F^"^(0)   =   [[[...[[F(u),A],  A]  ...  A]  =  ([  ,A]  )^'B 


From  (A. 2. 2)  and  (A. 2. 4)  we  get 
(A.2.5)      e-'^Be*  =  e"  '*' '  B 


Now  consider  the  operator 
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„r    N      -uA  u(A+B) 
G{u)  =  e    e  ^ 


(A. 2. 6) 

G(0)  =   1 


-uA   u(A+B)     -uA       u(A+B) 
G'  (u)  =  -  Ae    e       +  e    (A+B)  e 
(a  .2 .7) 

G'(u)  ^   e-^^  B  e^^^+B^   =   (e^^^  '^^^    B)    G(u) 


One  should  observe  that  the  operators  of  the  form 
gUU  ,    i  i-Q   fQjp  different  values  of  u  do  not  commute  with 
one  another. 

To  continue: 


u([  ,A]) 
(A. 2. 8)      G(u+du)  =  (1  +  du  e         B)  G(u)  . 


So: 


^   du  e^(f  >^^^B 
(A. 2. 9)      G(u)  =   TT   e 


For  our  purposes,  we  make  the  approximation  that 
[[[B,A],A],A]  f^  0. 
Then  we  get 


u([  ,A])  2 

(A. 2. 10)      e        B  ~  B  +  u[B,A]  + -^  [[B,A],A] 


Now  if  we  also  make  the  approximation  that 
B,  [B,A],  and  [[B,A],A]  commute  with  one  another,  then 
the  continuous  product  of  (A. 2. 9)  may  be  rewritten: 
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ou  2 

/       {B+u[B,A]+^   [[B,A],A])    du 

G(u)      %      e     °' 


(A. 2. 11) 

(uB   +^[B,A]    +^[[B,A],A]) 
G(u)     ^     e  ° 


Prom   (A. 2. 6)    and    (A. 2. 11)    setting  u  =  1,    we  get 

(A+B)      ^        A      (B   +-|Fb,A]    +-^[[B,A],A]) 


(A. 2. 12)  e 


~      e      e 
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Appendix  3 


00  4-ie 

(A. 3.1)      A^(x)  =  ^    r   exp  (  ^  +  Ixt)  dt 

-oo+le 

00  +le   CD  +ie 

(A. 3. 2)  A, (x)A,(y)  =^    T       T  dt  ds 

-oo+le  -oo+le 

3     3 
•  exp  (^  +  ^  +  ixt  +  lyt ) 


Substitute 


R   =  ^^^  u=t+s  t=|  +  v 


t- 3  c  "  „ 

r  =  x-y  V  =  -g —  3  =  2  -  V 


Then; 


00  +le 


(A. 3. 3)      A^(R+  |)A^(R-  f )    =  -^      /      du   exp    (  ^  +  IRu) 


-00  +le 

00 


/      dv  exp    (luv     +  Irv) 


-00 
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,2m 
(A. 3. 4)      ^o;^  A,(R+  ^)A,  (R-  ^) 


dr 


2m      1^'^      2'    1^''      2' 


00+le 


r=0 


m 


du 


-c»  +le 


-,  00 

•    exp    {^-  +  xRu)        v^"^  exp  luv^  dv 


but 


CD 


00 


(A. 3. 5)       /   dv  v^"^  exp    luv^      =     2      T 


o       /    ^        2m  .2 

2      /    dv  V        exp    luv 


-00 


-00 


Let  00  =:   e        '      uv     and   then 


00 


(A. 3. 6)     j^dv  v^  exp   luv^     -    ^^     ^^-^yg         y^ 
-00  u  0 


°°  m-1/2 


(D 


-O) 

e        doD 


00 


(A. 3-7)      /"dv  v^"^  exp   iuv' 


-00 


iir  ^  ,    1 
(e^  )        ^     „,    ^  l^ 
u        ' 


Thus 


,2m 


(A. 3. 8)      -^  A    (R+|)A    (R-  |) 
dr 


(-1)    e 


•(m+-i) 


r(m+  p 


r=0 


fir 


OD+le  .    3 

P         du   exp(-j^  +  IRu) 


u 


m+1/2 


-00+ie 


Let  u  =   ao); 
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,2m 
(A. 3. 9)      -^  A    (R+|)    A(R-  J) 


r=0 


_  ^(m+  ^)  ,        oo+le  ,    3  3 


,,    2     m-l/2 
HTT     a 


u 


12 

m+1/2 


-00  +le 


To  get  a   closed    from   for    (I. 5. 17)   and    (I.5.I8) 
of  the  main   text  when  D  is   an   odd   integer  we   set: 


(A. 3. 10) 


or 


3    ^"iwr 

^      "  2i~   ' 


R   =   "-,  y^       where         a   = 


aR  =   -    [E  -V(r)] 
2m[E-  V(r)] 


(mfi  VV) 


27T 


then 


(A. 3. 11) 


Itt/        In        +00 +le 

e  /'  dcD 

m+1/2   ^ 


lja3Vi„„i2 


2Trl 


/ 


l^tE-V(r))+i|52-|vv| 


-oo+le 
2m-l 


O) 


,m+l, 


(-l)'"'^^27r(2mh|VV|)     ^  ^2m 


''^^ 

(2m)~^~     r(m  +  h 


dr 


2m     l^p73       2^^i^p75       2^ 


r=0 


D+1 


Thus   using    (I.5.I7)    of  the  main  text  we   let  m  =  — p—  and  get 


D+1 


(A. 3. 12)    p{r) 


{- 

-1)' 

2 

hir  Ojj(2mR 

VV  )^/ 

5        ,D+1 

Dh'' 

du^-^1 

^l\l/5 

f  -)    A 

\l/3   -   2^ 

D   =    1 

,5,5,    . 

•    •         • 
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u=0 


And  using  (I.5.I8)  of  the  main  text  we  let 
m  =  — p-^  and  get: 


(A. 3. 13)  t(r)  = 


D+5 
(-1)    ^      kir  0^(2mh 

D+2 
5 

^D+3 
du^+5 

(D+2)    2in   h^ 

^i^l 

-a 
,1/3  - 

2' 

u=0 


D  =  1,3,5,  .-.  . 

While  we  are  at  it  we  notice  that  by  comparing  (I. 5. 28) 
of  the  main  text  with  (A. 3. 12)  we  see  that: 

2k-l 


a      '^K-l 
(A. 3.1^)  r  (a  -t)  ^   A^(-t)  dt 


-00 


2k-l 


=  (-1)^  2v   2 


,2k 


du 


5Kf*i(W3-^>'^i'7i75 


|)] 


u=0 


for  k  =  0,1,2, .. .  . 
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Appendix  4 

Asymptotic  Expansions  of  the  mass  density  p(r)  and 
kinetic  energy  density  t(r). 

From  the  main  text  let  us  rewrite  equations  (I. 5. 17) 
and  (1.5.18) 

lirD    00  +le 


(I. 5. 17)  p(r)  =  -  ^  (2m)^/2  r(|)  I 


h' 


^       r     dco 

D  '"""  '  ^2'    27rl      /     D+2 


-00+ie 


lao^h^  i„.ri2 


-lco(E-V(r))+^|^  |VV| 
•  e 

1)  ,^  ^  2   ^,D+2^  e         r    dcD 


(1. 5. 18)  t(r)  =  -  -^  (2m)  "   r(^)  ^ ^ /   ^^ 


h  27rl     -     -2- 

-00  +ie 

-l^(E-V(r))-f^  1VV|2 


CASE  I  (Oscillatory  Region  E  >  V(r)); 
Here  we  substitute 

3      24m   '~ 


u6   =  cd(E-  V(r)) 

and  get  ^^^ 

r(D/2)ru   2m(E-V(r))   D/2   ^""T" 


+00  +le  -, 

J    -^   exp  [1^(^  -  u)] 


-oo+i  e  ^ 
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(A. 4. 2)      t(r)    = ^-^   [   ^-^^   ]    2        ^ 


2m  h 


P  P  ^  27ri 

00  +le  ^ 

.y'      ^g^     exp    [ip(-^   -   u)] 

-00  +ie 


where 


3/2 
[2m(E  -V(r))] 

(A. 4. 5)  P      =     - 


mh    |VV| 

It   is    clear   that   to   find   the  asymptotic   expansions 
for    (A. 4.1)    and    (A. 4. 2)    we  need   only   consider: 

IttD        00  +le 
(A. 4. 4)  Jd  =   -  fFi"         J        -^     exp[ip(^-u)]     . 

-oo+le  "^ 
To  get  the  path  of  steepest  descent,  let 

h  (u)   =  1  ( J-  -  ^) 

h'(u)   =  1  (u^  -  1) 
So  we  get  the  saddle  points  u  =  +  1. 

Letting  u  =  X  +  ly  we  find  that  the  equation  I^(h(u)  )=Ij^(h(l) ) 
corresponds  to  the  equation: 

2     (x-l)^(x+2) 

y    =  5x — 

and  the  equation  I  (h(u))  =  I  (h(-l))  corresponds  to 


the  equation: 


2     (x+l)^(x-2) 

y    =  — 5r — 
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We  thus  get  the  following  graph  for  lines  of  steepest 
descent: 


TT 


Here  the  arrows  point  in  the  direction  of  steepest  descent. 
Thus  we  choose  for  the  path  of  integration: 
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u  =  X  +ly 


We  now  designate  the  contribution  to  Jj^  from  the 
path  segment  P,  as  J^. 

Now  Jj^,  Jj^,    J^   and  J-^  may  be  evaluated  by  standard 
methods  [I3] . 

Carrying  out  the  computations,  the  author  got  the 
following  answer. 
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Jd  +  Jd  +  Jd  +  Jd    - 


sln(f  .  ^d^) 


TT 


00 


n 


(-1)^'  r(6n  +^) 

n-n  U^^^      2n+l/2 

"-°    (4n):    3  "^  P 


> 


^_       -        (-l)"r(2n.|)    Ag^^^ 
D+2  ^      ^— 7T  ^2n+  5/2 


(A. 4. 5) 


r(^)    ^^ 


P 


+ 


cos(#  +  :^^^)   r  ^      (-1)"  r(6n  +  i) 


3 


TT 


> 


L"=0    (4n4-2)I    3^n+2p2n+3/2 


+ 


rD+2. 


z 


^      (-1)"  r(2n+|)    A^^^^ 


r(^)    7^ 


p 


2n+  5/2 


where 


,D 


n 


n  4 — 


n 


r( 


3n-3k+l 


) 


^=0     3"-^(n-k)l    r(^^^4^)(k+l)! 


Z 


k    r(k-^+§  +2)  r(i+  ^  +^) 


2    ,^/    X  vx,  ,    2    ■    2/    ,    k 


£=0 


!>'  r(^) 


For  J^,    J-p,    J^     set   e   =   l/p 


For  J^^,    let  u  =  e^^^^^^o)  =  e^^^^^  f  and 
;5      1^+-   ,,   _   ^-l7r/2   _,^  _   ^-lTr/2  (^ 


for  JfJ,    let  u  =  e   ""'  "   ecD  =   e   ^"'  "  —   .      Then   we  get 


P 
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^D/2      .      ttD         OD  ^2 

-z  I-  S  Sm   -?r-  p  -CD  36 

(A.^.6)     4-^1    -     - ^-^   J      -     ^l,  . 

1  od"^ 

„         T.4    ,     ,  10         1      10 

For  J^letu  =   £e        =-^e 

130 
D/2     ^     ^2  -  ^^      -le^®      ^   ^   2 

(A. 4. 7)      Jp     =     ^— 2^j:^ f    6Q   e        ^      e  e      ^ 

-7r/2 


3  S 

When  D   Is   even,    J^  +  J^  =  0  and 

J,         _D/2    ^"~ir  p         ^  OD      ,    ,^k      Ik© 

(A.H.8)      4=^ ^—        /      e      2      (^lzi)_e  ) 

CD         ,..£  13-^© 

^=0  i:(3p  ) 


and  thus  for  D  even   noting  that  the  only  terms  that 
contribute  to  J^  are  those  for  which  k  =  D/2  -  3^, 
we  get 

i^    D/2  [D^] 
(A. 4. 9)   Jt^  +  J^  +  J^  =  P     XZ  -n "^ pT  ^^  ®v®"^ 


For  D  odd  things  become  somewhat  more  complicated. 
From  (A. 4. 7)  we  get 

,  2p^/2   sin  ^         00  ,    .a 

(A. 4. 10)      J^  -      -     ^  ^~        ^"-^A   p 

^  (-l)^__    _ 


k=0     k![6i+2k-D] 
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3     R  2 

To  compute  Ji.   +  J^  let  us  substitute  w  =  u  and  we  get 

OoD/2     .     TtD    (D  rs 

IT  ^ 

^   (-1)^  U^^   ^ 


Now  if  we  make  use  of  the  equation 

°°   .2        _rz  ..2  _ra_     ^k  2k+l     ^ 


(A. 4. 12)   Z'e"^   dt   =  ^  -  e"^' 


it  is  possible  to  show  that 


r  2N  -t^       r(N+|) 
(A. 4. 15)   /  t"^^  e  ^  dt  =  ^-^ 


n=0  1.3-  •••(2k+l) 


2  00       ok  2N+H-2k 

_  e-^  2r  — ^-^ 


k=0  (2N+l)(2N+3) . , .(2N+l+2k) 
thus 

(A.4.14)   J^  +  J5 _ 2_ 

f-        (-l)\(p(3,_  D)_2,-l  ^  2?^ ) 

ra  ^!(3P^)        ^       k=0  (6i-D)(6i-D+2)...(6^-D+2k)  _ 


127 

Handbook  of  Math.  Functions,  edited  by  Milton  Abramowltz 
and  Irene  A.  Stegun  (New  York,  Dover  Publications,  Inc.  I965), 
p.  297  (7.1.6). 
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On  the  other  hand  one  can  show  by  multiplication 
of  series  that: 


(A. 4. 15)   e^   ^  -L^)   ^ 


k=0   k;(2k-D+6^) 

fco         'H^   U-n)  •n:{2k-D+6^) 
And  by  Induction  on  k  we  have: 

(A  h  16)   y^  (-1)  ^ 2 

^  Ck-n)In!(2k-D4^i)   T6i-D) {6i-D+2) . . . (6i-D+2k) 
Setting  X  =  1  it  is  now  clear  that: 


-1  ,0°  o^ 

(A. 4. 17)   e    ZZ  (bii-ij)(bi^-D+2)...(6i-D+2k) 


=  ? 


^      f-1)^ 


^  kI{b^+2k-D) 

Now  using  (A. 4. 17)  to  add  (A. 4. 10)  and  (A.4.14)  we  arrive 
at  the  expansion: 

„D/2   .   ttD  - 

(D  odd) . 
Finally  using  the  relation 
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r(5^-  |) 


IT 


ly 


r(|  -3i+l)    sin  (3i-  |)Tr 


we  now  see  it  Is  possible  to  combine  (A. 4. 9)  and  (A.4.18) 
and  arrive  at: 


(A.i^.19)  4+ 4  + jI-^^^^'z 


00 


-D 


2si 


^=0  r(^  -5i+i)ii(3P  ) 


(D  an  arbitrary  Integer) 


Thus    from   (A. 4. 5)   and    (A. 4.19) 


,D/2 


00 


1 


'^'"^         iS  r(|  -3n-H)nI(3p^)" 


(A. 4. 20) 


TV 


1 


^    (-l)'^r(6n+^)  ^  CO     (-l)"r(2n+|)Ag^_^^ 


fco    (4n)I3'^¥"-*-^2  ■  r(^)  ^  ,2n+5/2 


3 


4-cos(^+^li^) 


TT 


(-lfr(6n+|)  ^  00     (^fV(2n+|)A^ 

(4n+2)l3'"-^2p2n4-3/2   +  :^^  gl  2HW2 


P 


where 


,D 


^     =Z 


n 


r(^^l=fS±i) 


n 


and 


^=°  5"-k(n-k).T(^i^)(k+l)l 


k    r(k-£-f^+2)r(^+|+^)    , 


Z 


P  = 


[2m|E-V(r))]^/2 
mft    |VVI 


ly 


Ibid. .    p.    256    (6.1.17) 


98 


Prom  (AA.I)   and  (AAA) 


(A. 4. 21)    p(r) 


2m(E-V(r) ) 


D/2 


Jj5   (E  >  V(r)) 


and  from  (A. 4. 2)  and  (A. 4. 4)  we  get 


(A. 4. 22)   t(r) 


r(^±^)o 


2m  h 


D 


2m(E-V(r) ) 


D+2 


J. 


D+2 


Case  II   (Exponentially  Decreasing  Region  E  <  V(r)): 

Here  we  modify  (I. 5. 17)  and  (I.5.I8)  of  the  main 
text  by  substituting: 


3  24~m 


up'      =     a)(V-E) 


thus: 


(A. 4. 23)    p(r)    = 


-r(|)o^ 


h 


D 


2m(V(r)-E) 


D/2 


IttD 


+oo+ie 

r  du 

27ri         J  D+2 

-00  +le     u  '^ 


exp  i  6  '  (^  +  u ) 


and 
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(A. 4. 24)      t(r) 


I  \    2    ^     D 


where 


D 


2m  h 


D+2 


2m(V(r)-E) 


D+2 


(e"") 


+00  +le 


du 


27rl  J  D+4 

-00  +ie  2 

u 


exp    (ip'(^  +u)) 


[2m(V(r)-E)]^/2 


P'    = 


mtr   vv 


Here  we  proceed   similarly  to   Case   I  and    consider   only 
the   integral 

l-irD 
(A.4.25)      I 


+00  +ie 


,3 


'D 


Now 


2Tri         J  D+2 

-00+ie  2 


h(u)  =  l(-:=—  +  u) 
0 


^"^        exp  13 '(y-  +  u) 


h'(u)  =  i(u  +  1)  . 

So  we  get  the  saddle  points  u  =  +  i.   Letting  u  =  x  +iy 
we  find  that  Im(h(u))  =  Im(h(+1))  =  0  corresponds  to 
the  equation: 

x(x^-  5y^  +5)  =  0  . 

Prom  this  we  get  the  following  graph  for  the  paths  of 
steepest  descent. 
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\ 


> 


(0,-1) 


Here  as  in  Case  I  the  arrows  point  In  the  direction 
of  steepest  descent.   Now  we  see  that  our  path  of 
integration  must  be: 


y  f 


0  = 


TT 


-> 


X 
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The  rest  of  the  computation  is  carried  out  by 
standard  methods  [13]'   My  result  was: 


{AA.26)      Ip  '^' 


-23' 

2tt 


00   (-1)  r( — 2~' 

n=0  (2n)  13  p'^   '^ 


r(5±2)  S=o 


f-  ^(^'t-i'"*L.i 


,n+3/2 


whereas  in  Case  I 


n 


P' 


n 


r(^nz|k±l) 


k  r(^^+i5+Mr(k-i+2+g) 


5  (n-k).'3"-^(^^)  T- 
[2m(V(r)-E)]5/^ 


:0 


2 

(k+1)]  r(^)  3 


mK  |VV| 


From  (A. 4. 23)  and  {kA.2^) 


(A.i^.27)   p(r)  = 


2m(V(r)-E) 


D/2 


■I 


D 


and  from  (A. 4. 24)  and  (A.il.25) 


D+2 


(A. 4. 28)   t(r)  = 


i  ^  2  ^  D 


2m  h' 


D 


2m(V(r)-E) 
P' 


-D+2 
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